1 



1. Introduction 

It is well known that the far infrared fine-structure lines of 
abundant elements (Oxygen, Carbon, Nitrogen, Silicon, Sulfur, 
etc), either in their neutral or ionized states, contribute most 
of the gas cooling of the interstellar medium in galaxies 
dHollenbach & Tielens Il999t iGoicoechea et al.ll2005h . Far in- 
frared fine-structure lines from ionized atoms are useful trac- 



coo Ung is at least similar to the contributio n of the ground state 
line dStutzki et al.lll997tlBavet et al.ll2004 . 
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ers of HII regions (e.g. IGoicoechea et al.l 120051) . For the neu- 
tral ISM, the main cooling lines are those of ionized carbon 
[CII] and neutral oxygen [OI]. By a nalyzing a set of ISO obser- 
vations of external galaxies, iMalhotra et al.l (120011) have con- 
cluded that PDRs contribute to a large fraction of the emis- 
sion of the [CII] and [OI] lines. In molecular gas, the cool- 
ing radiation due to atomic carbon and carbon monoxide is 
significant. The theoretical predic tions dGoldsmith & Langei 
ll978HHollenbach & Tielenslll999h have been confirmed by the 
COBE-FIRAS and ISO-LWS observations of the Milky Way : 
apart from [CII], [OI], [Oil] and [Nil], the most intense Unes 
are from C and CO. The relative contributions of the different 
lines of C and CO vary along the Galactic plane ; the bright- 
est lines and more excited states are seen towards the Galactic 
Center, while the rest of the disk shows lower excitation. Also, 
C contributes less in proportion to the t otal cooling towards 
the Galactic Cente r than towards the disk dBennett et al.lll994 ; 
iFixsen et al.lll999l) . 

The CO cooling is typically provided by the submillime- 
ter Unes, with rotati onal quantum numbers between 3 and 
8 dBavet et alj 120041) . Such observations can be performed 
from high altitude dry mountain sites, such as the Mauna 
Kea summit in Hawaii. Previous studies using the James 
Clerk Maxwell Telescope (JCMT) and the CSO have shown 
that the '~CO(4-3) l i ne is generally br i ght in galaxy nuclei 
dGuesten et al.1 1 19931; llsrael et al.l [l99l llsrael & Ba3 12001 
120031: iKrameret al.lbOOsl ...). Much less information is avail 



able on the other high frequency CO lines. The first re- 
port of extragalac tic '^CO(6-5) detections was published by 



Harris et al.l(ll991b 15 years ago. In the past years new analyses 



of the '^CO(6-5) and '^CO( 7-6) emission have been presented 



by IWard et al.1 (120011 120031) f or M 82, and by [Bradford et al 
d2003h and lBayet et al.l (l2004l) for NGC 253 and Henize 2-10 



A first co mparison of the CO line spectrum, also called "CO 
SED" (see IWeiB et al.ll2005l: iMao et al.jl200d ...) in star form- 
ing galaxies is presented bv lBavet et al. d2004l) . They show that 
the CO line spectrum is very similar in the two star forming 
galaxy nuclei, NGC 253 and Henize 2-10. Distan t starburst 



galaxies seem also to share the s ame CO spectrum dCox et al . 
2002; Bertoldi et al. 2003; CarilH et al. ll20()4l: IPetv et al.l l2004l: 
Walter et aLkO^ICarilh et al.ll2005 



...). 



Although contributing less than CO, atomic carbo n is an 
important coolant in the ISM dOerin & Phillipsll2000l) . While 
observations of the ground state line (^Pi - ^Po at 492 GHz) 
can be found in the literat ure for a few te ns of sources (e.g. 



Gerin & Phillipsll2000; Israel & Baasll2002l) there are very few 
reported detections of the excited line (^Pa - ^Pi at 809 GHz) 
although the contribution of the latter line to the neutral carbon 



As potential tracers of the gas cooling, submillimeter C and 
CO lines are expected to provide information on the gas heating 
rate, which is dominated by the incident FUV radiation, mainly 
due to massive and young stars. Therefore, the molecular cool- 
ing lines are expected to provide information on the galaxy star 
forming activity, as do the fine structure lines in the far infrared. 
In order to have a full picture of the CO cooling, the contribu- 
tion of missing CO lines (blocked by the Earth's atmosphere) 
can be predicted from the series of observed lines, using state- 
of-the-art radiative transfer models. The method has been pre- 
sented in our previous paper (iBayet et al.ii2004i) . In this paper, 
we present in this paper results of a survey of the C and CO sub- 
millimeter lines in a sample of nearby galaxies. The data are 
used for the two following purposes : /) determination of the 
molecular cooling rate in galaxies of different morphological 
type; ii) study the shape of the CO cooling curves obtained in 
the target galaxies. From this analysis, we show that the com- 
bined information on C and CO submillimeter hne spectra can 
be used as a powerful diagnostic of galaxy star forming activity. 

The galaxy sample is presented in section 2, the observa- 
tions parameters are described in section 3, while the resulting 
spectra and maps are introduced in section 4. We discuss in sec- 
tion 5 how we use LVG and PDR models for fitting the series of 
observed CO lines. In section 6, we compare results obtained 
for the center of Milky Way and for the Cloverleaf QSO with 
those derived from this work for our galaxy sample. The main 
conclusions are summarized in section 7. 



2. The sample 

We have selected galaxies which are bright in the '^CO(l- 
0) and '^CO(2-l) lines and are nearby (distance less than 20 
Mpc, except for the two ULIRGs Arp220 and Markarian 231). 
Galaxies have also been selected to have a large variety of 
galaxy types. The sample includes normal spiral galaxies (IC 
342, M 51, NGC 4736, NGC 6946), starburst galaxies (M 82, 
M 83, NGC 253, NGC 3079), irregular, star forming galax- 
ies (IC 10, Henize 2-10), interacting galaxies (The Antennae, 
NGC 6090), ULIRGs (Arp 220, IRAS 10565+2448, Markarian 
231) and the elliptical galaxy Centaurus A. In this work, we 
analyzed two positions in the Antennae galaxies : the nucleus 
of the northern component, NGC 4038, and a position named 
"Overlap" hereafter, which corresponds to the position of the 
most massive H2 concentration, not very far from the nucleus 
of NGC 4039 (The coordinates of NGC 4038 and Overiap are 
listed in Table [U. Intense MIR emission due to star formation 
has been detected at the Overlap positi on ( Vigroux et al ll996l) 
as well as bright CO Unes (GMC4-5 in lWilson et al.ll2000h . 

Properties of the sample galaxies are summarized in Table 
[1] Although it is not a complete sample (because these obser- 
vations are difficult and time consuming) it includes represen- 
tative types of nearby galaxies. 
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Table 1. Basic properties of the sample galaxies. 







RA('1950~) 


DECn950') 


Dist. 
(Mpc) 


Velocity 
LSR 

(kms-') 


Optical 
size" 


Metallicityi 


IC 10 


dirr IV/BCD" 


00:17:44.0 


59:00:18.0 


1' 


-344 


6.8'x5.8' 


8.31+0.20'* 


NGC 253'' 


SAB(s)c;Hii, Sbrst 


00:45:05.7 


-25:33:38.0 


2.5'** 


240 


27.5' X 6.8' 


8.99± 0.31"' 


IC 342 


SAB(rs)cd HII" 


03:41:57.2 


67:56:27.0 


1.8^ 


35 


21.4'x20.9' 


a9.30''' 


Henize 2-10'' 


10 pec, Sbrst 


08:34:07.2 


-14:26:06.0 


6''' 


850 


30" X 40 " 


-8.93"' 


M 82 


I0;Sbrst HII" 


09:51:43.8 


69:55:00.9 


3.2^ 


200 


11.2'x4.3' 


9.00+0.12'* 


NGC 3079 


SB(s)c;LINER Sy2" 


09:58:35.0 


55:55:15.4 


15.6* 


1331 


7.9'xl.4' 




IRAS 10565+2448 


LINER HII" 


10:56:36.2 


24:48:40.0 


172^ 


12923 


0.4'x0.3' 




NGC 4038 


SB(s)m pec" 


11:59:19.0 


-18:35:23.0 


13.8^ 


1634 


5.9'x3.2' 




Overlap 


SA(s)m pec" 


11:59:21.1 


-18:36:17.0 


13.8^ 


1510 


3.rxl.6' 




NGC 4736 


(R)SA(r)ab;Sy2 LINER" 


12:48:32.4 


41:23:28.0 


4.3^ 


314 


11.2'x9.1' 


9.01+0.17'* 


Mrk231 


SA(rs)c?pec Syl" 


12:54:05.0 


57:08:39.0 


173.9** 


12650 


1.3'xl.O' 




Centaurus A 


SO pec Sy2" 


13:22:31.6 


-42:45:32.0 


3.5'^ 


550 


25.7'x20.0' 




M51 


SA(s)bc pec; HII Sy2.5" 


13:27:46.1 


47:27:14.0 


9.6'° 


470 


11.2'x6.9' 


9.23+0.12'* 


M83 


SAB(s)c; HII Sbrst" 


13:34:11.3 


-29:36:42.6 


3.5" 


516 


12.9'xll.5' 


9.16+0.12'* 


Arp220 


S?;LINER;HII Sy2" 


15:32:46.7 


23:40:08.0 


77'" 


5450 


1.5'xl.2' 




NGC 6090 


Sd pec HII" 


16:10:23.9 


52:35:11.0 


118'^ 


8831 


2.8'xl.5' 




NGC 6946 


SAB(rs)cd HII" 


20:33:48.8 


59:58:50.0 


5.5'^* 


50 


11.5'x9.8' 


9.06+0.17'* 



References: 1: Adopted value (se e text) ; 2: iMcCalll h989l); iKarachentsev & TikhorioyI (Il993l) ; 3: iDumke et alj feOOlh consistent with 
the value of iTammann & Sandage'. l ll968h; 4: ISofue et al.l lll999l); 5: lOlenn & Hunted (200 ih; 6: ISaviane et al " 



(119871); 8: 



. iBrvant & Scovillei 11999i) ; 9: Ide Vaucouleurs et al 

iGerin & Phillips! (Il999h with Hp = 7 5 kms 'Mpc '; I3 :|Tulb! 
( Il998l) ; 16:|Zaritskv et 31.71994; 17: |johanssoij i fl987h ; 18 



19791) ; 10: Ide Vaucouleurs! (Il99lh ; 1 1 : iThim et al 




7: ITuUv & FisheJ 
12: Redshift from 



\; 14:lArimoto et al.l (ll99^15:IVila-Costas & Edmunds! (ll992l) : lGamet3 



database; : See iBavet et al. C2004 



As in lMauersberger et alj ( 1 19961) (see lBavet et al.l ( |2004|) ); 
to obtain more information on the properties of NGC 253 and Henize 2-10. 



Data from the NED 



3. Observations 

The observations were made during various sessions at the 
Caltech Submillimeter Observatory (CSO) in Hawaii (USA) 
with the Superconducting Tunnel Junction receivers operated 
in double-side band mode. The atmospheric conditions varied 
from good (T225 S 0.1) to excellent (T225 ~ 0.06). We used a 
chopping secondary mirror with a frequency of around 1 Hz. 
We used a 3' chopping throw for the [CI](¥i-¥o), CO(2-l), 
CO(3-2) and CO(4-3) lines. There is no sign of contamination 
by emission in the off beams. We restricted the chopping throw 
to 1' for [Cl](^P2-^Pi), CO(6-5) and CO(7-6) as the emission 
is very compact in these lines. Spectra were measured with 
two acousto-optic spectrometers (effective bandwidth of 1000 
MHz and 500 MHz). The first one has a spectral resolution 
about 1.5 MHz and the second one about 2 MHz. The IF fre- 
quency of the CSO receivers is 1.5 GHz. The main beam effi- 
ciencies (77) of the CSO were 69.8%, 74.6%, 51.5%, 28% and 
28% at 230, 345, 460, 691 and 806 GHz respectively, as mea- 
sured on planets'. For the [Cl](^Pi-^Po) and the [Cl](^P2-^Pi) 
lines, we used receivers at 492 and 809 GHz, so = 51.5% and 
28%, respectively. We used the ratio - to convert T^ into T„,h. 
The beam size at 230, 345, 460, 69 f and 806 GHz is 30.5", 
21.9", 14.5", 10.6" and 8.95"Q. The pointing was checked us- 
ing planets (Jupiter, Mars and Saturn) and evolved stars (e.g. 



IRC 10216, R-Hya, CRL 2688, CRL 618, NGC 7027, R-CAS 
and O-Ceti) for all lines except CO(7-6). Planets were the sole 
pointing sources at 806 GHz. The pointing accuracy is around 
5". The overall calibration accuracy is ^ 20%. Data have been 
reduced using the GILDAS/CLASS data analysis package. The 
spectra have been smoothed to a velocity resolution of 10 
kms Gaussian profiles have been fitted to observed spectra 
(see Figs. [T]to [T6T l, and linear baselines have been subtracted. 



4. Data analysis 

4.1. Spectra and maps 

Spectra of the galaxy nuclei are shown in Figs. [T]to [141 The 
spectra taken at positions in the spiral arms of M 83 are shown 
in Fig. [15] while those taken at positions in the NGC 6946 arms 
are presented in Fig. [161 The observed line is indicated above. 
Figs. [TSland [TTlshow integrated intensity maps of the '^CO(3- 
2) (top) and the '^CO(6-5) (bottom) lines towards M 83 and IC 
342, respectively. 

Tables IB. II and IB. 21 in appendix [Bllist the line intensi- 
ties (A in Kkms"' and 1 in Wm^^sr"') and the line fluxes (F 
in Wm"^) resulting from Gaussian fits for the sample galaxies 
(for each observation, we give the corresponding b eam size) . 1 



(in Wm ^sr ') is derived using Eqs. (3) and (4) in lBavet et al 



See web site: |http://www.submm.caltech.edu/cso/| 



(120041) . To compute the flux, F (in Wm ^), we used Eq. (5) 
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IC 10 



(0.0) 




^ '00(6-5) 



CK^Pj-^P,) & '^C0(7-6) X 5 



(0,30) 




-400 



-350 
Velocity (krm/s) 



-300 



Fig. 1. Observed spectra towards IC 10 nucleus (see Table IB. II ) 
except for the [CI](^P2-^Pi) and the '-CO(7-6) spectra which 
correspond to the offset position (0",30")- Velocities (horizon- 
tal axis) are given in kms"' relative to the LSR (Vlsr) and the 
line intensities (vertical axis) are in units of T,,,^ (K). The grey 
curves are Gaussian fits. The observed line is written above 
each spectrum. In this figure, we present the [CI](^P2-^Pi) and 
the '^CO(7-6) lines on the same spectrum since they have been 
observed simultaneousl y. The [CI](-'Pi--'Po) spectrum is from 



Gerin & PhilUpsI ( |2000|) but it has been analyzed again to ob- 
tain an homogeneous dataset of observations. 



Bavet et al.l ( 120041) . The estimated errors on data listed in 



Tables IB. H and IB. 21 are also indicated in column 5. For most 
sources, t he atomic carbon [Cl1( -^P|--^Pn) data have been pub- 
lished by iGerin & PhillipsluOOd but they have been analyzed 
again for consistency. All the CSO observations we obtained 
for our galaxy sample are summarized in Table IB. II where we 
list the data for all positions in the nucleus as well as in the 
spiral arms when it is appropriate. 




-200 



-100 



Velocity (km/s) 



100 



200 



Fig. 2. Observed spectra towards IC 342 nucleus (see Table 
IrTT ). See the caption of Fig. □ We present the [CI](¥2-^Pi) 
and the '^CO(7-6) lines on the same spectrum since they have 
been observed simultaneously. 




100 200 
Velocity (km/s) 



Fig. 3. Observed spectra towards M 82 nucleus (see Table IB. It . 
See the caption of Fig. [T] 



4 



1 L_l' 



NGC 3079 ■ 





1000 1200 1400 

Velocity (km/s) 



Fig. 4. Observed spectra towards NGC 3079 nucleus (see Table 
IB. It . See the caption of Fig. [1] 




100 

Velocity (km/s) 



Fig. 5. Observed spectra towards IRAS 10565+2448 nucleus 
(see Table IB. 11 1. See the caption of Fig. [T] 



Table |B.2| in the electronic appendix presents central po- 
sition data for a restricted sample (IC 10, IC 342, M 83, NGC 
4038, Overlap and NGC 6946) we have been able to observe 
extensively (up to the '^CO(6-5) or '^CO(7-6) line). For these 
sources, we list the available informations in the literature. 



in the same way as iBayet et al.l (l2004h did for NGC 253 and 
Henize 2-10. All CO transitions we used for the modelling 
work, are identified with an asterisk in column 8 of Table IB. 21 
Moreover, in order to compare at the same spatial resolu- 
tion different CO line intensities, we have convolved the line 
intensities to a common (final) beam size of 21.9", which is 
the CSO beam size at the frequency of the '^CO(3-2) line. 




X0(7-6) X 2 



—I I I l_ 



1500 



1600 1700 
Velocity (km/s) 



1800 



Fig. 6. Observed spectra towards NGC 4038 nucleus (see Table 
IB. II ). See the caption of Fig. [1] 



Precisely, to perform this convolution, we multiplied A, I and 
F by factors depending on the size of each emitting source (See 
Eqs. lA.ll to IA.4l in appendix |A|, the initial (observed) beam 
size and the final beam size (21.9"). To determine the source 
sizes, we used high spatial resolution maps and fitted Gaussian 
spatial profile to these maps. The sources may be either ax- 
isymmetric (a x a) or elliptical (a x b). Table |2]lists the adopted 
sou rce sizes and the m aps we used for this measurement. 



Bayet et al.l (l2004l) used the same method in their study of 



NGC 253 and Henize 2-10. 



4.2. The C and CO cooling rates 

We have derived the observed C and CO cooHng rates which 
provide an essential information for the thermal balance of the 
studied galaxy nuclei. We estimated the observed C and CO 
cooling rates by summing the intensities (I in Wm^^sr"') of all 



CO transitions listed in Table IB. 21 and identified with an as- 
terisk (both literature data and our dataset). We have computed 
the observed C and CO cooling rates in the galaxy nuclei for a 
common beam size of 21.9". This corresponds to linear scales 
of 106 pc, 191 pc, 1.5 kpc, 372 pc and 584 pc for the IC 10, IC 
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Table 2. Source sizes used to convolve the data to 21.9". For each source, we used maps (specified in column 5) to determine the 
shape of the Gaussian intensity profile (axisymmetric : a x a or elliptical : a x b) and estimate the values of the full widths at half 
maximum (FWHM) 





FWHM 


line (frequency) 


Observed beam size 


References 


IC 10 


20.0" X 20.0" 


'-CO(4-3) (461 GHz) 


14.5" 


This work' 


NGC 253" 


23.0" X 11.0" 


'2CO(6-5) (691 GHz) 


10.6" 


Bavet et al. (2004)- 


IC 342 


11.0" X 11.0" 


'^CO(6-5) (691 GHz) 


10.6" 


This work' 


Henize2-10'' 


13.0" X 13.0" 


'^CO(3-2) (345 GHz) 


22" 


Meier et al. (2001 f 


NGC 4038 


13.2" X 9.9" 


'2CO(1-0)(115 GHz) 


4.91" X 3.15" 


Wilson et al. (2000)' 


Overlap 


1 1.0" X 8.8" 


'^CO(1-0)(115 GHz) 


4.91" X 3.15" 


Wilson et al. (2000)^ 


M83 


12.0" X 10.6" 


'-CO(6-5) (691 GHz) 


10.6" 


This work^ 


NGC 6946 


17.0" X 17.0" 


'2CO(4-3) (461 GHz) 


17" 


Nietenetal. (1999)** 



We include the sources from iBavet et alj (120 04^. Reference s: :The 20" FWHM is in agreement with the size deduced from the CO maps 



found in Wilson & Re id (199jJ) or in iBolatto et al. ( 200(1); ^ : iBavet et al.l ( |2004|) found good agreement with the source size d educed from 



IPeng etalj ( ll996h map; ^ : The 11" FWHM is very similar to the size deduced from the '-CO(4-3) map of llsrael & BaasI |2003h.lEckart et alj 
1 19901) using the lRAM-30m telescope, obtained a map in the '^CO(2-l) line wh ich shows a source size around 20".; : See lBavet et alj ( l2004l) ; 
^ : OVRO map; ^ : Overlap corresponds to SGMC 4-5 in this OVRO map; ' : iLundgren et"ZI ||2004|) considered the M 83 nucleus as a 12" 
axisymmetric source which agrees with our value; * : The '-CO(l-O) map published bv lWeliachew et al] jl988l) . is consistent with a 15" source 



342, NGC 4038 (and the Overlap region), M 83 and NGC 6946 
nuclei, respectively. The observed cooling rates for C and CO 
are listed in Table [3]together with the CO cooling rates derived 
from PDR and LVG models (see Sect.lS^. 

For all studied galaxies, the CO lines contributing the most 
to the observed CO cooling rates are '-CO(6-5) and '-CO(7- 
6), followed by '2cO(4-3) and '^00(3-2), with varying relative 
contributions of those lines in the studied sources. These results 
confirm that the CO transitions with high-J {i upper > 3) are con- 
tributing the most to the total observed CO cooling rate, with 
the highest contribution for the '-CO(6-5) line in the observed 
galaxies. 

For all targets, it is also noticeable that the observed CO 
cooling rate is higher than the observed C cooling rate, by a 
factor of < 4.0 for IC 10 and IC 342, of 6.9 and 19.7 for NGC 



253 and Henize 2-10, respectively (see iBavet et alj|2004l for 
these last two sources), of >20 for NGC 4038, Overlap, M 83 
and NGC 6946, respectively. In the NGC 4038, Overlap, M 83 
and NGC 6946 nuclei, the line which contributes the most to 
the C cooling rate ([CI](^P2-^Pi), as it is shown in the follow- 
ing) has not been observed yet, but the difference between the 
C and CO cooling rates is large enough not to modify the dom- 
inant role of carbon monoxide for the gas cooling with respect 
to the atomic carbon. 

Similar results have been observed for the distant 
galax ies Jl 148+5251 (z^6.42) and PSS2322-H1944 (z=4.12) 
(See ICox et alj l2002t iBertoldi et al.1 [20031: IPety et al.l 12004 : 



Walter et al.l 2004 ). The lines which contribute the most to the 



CO cooling rates are those with Jupper^ 3. The difference be- 
tween the C and CO cooling rates seems also to be important 
(factor of > 10) in distant objects. 

Because the observed C and CO cooling rates are computed 
with the observed lines only (see asterisks in Table |B.2| i, unob- 



served Hnes (e.g. [CI](^P2-^Pi), ^2CO(5-4) or ^^CO{8-7)) might 
contribute significantly to the CO or C cooling rates. Therefore, 
we used models in Sect. |5]to predict the missing line intensi- 
ties; and also to obtain the physical properties of the warm and 
dense gas contained in the galaxy nuclei. Table [3] shows that 
the observed CO cooling rate is x 70% of the modelled CO 
cooling rate for galaxies with '^CO(7-6) detections and 25%- 
50% for galaxies without '^CO(7-6) detections. 

5. CO models 

5.1. Description of tine models 

In this section, we used the measured CO line ratios from the 
line intensities (I and A) to determine the physical conditions 
of warm and dense molecular gas in galactic nuclei; namely 
the kinetic temperature (Tk), the gas density (n(H2)), the CO 
column density divided by the line width (N(CO)/Av) and the 
Far UV radiation field (xfuv)- In the first section, we use a LVG 
radiative transfer model and in Sect. 15.31 we discuss the use of 
a PDR model. XDR model results applied to AGN nuclei from 
our sample will be presented in a forthcoming paper 

The radiative transfer model s, based on the LVG for - 
malism, have been developed by Goldreich & Kwaiil ( 1974 ): 
de Jong et alJ (119751) . The source is modelled as a one com- 
ponent spherical cloud, with uniform kinetic temperature and 
density. When using both '^CO and '^CO data, there are four 
main variables in LVG models : the molecular hydrogen den- 
sity n(H2), the gas kinetic temperature Tk, the CO column 
density divided by the line width, and the '^CO/'^'CO abun- 
dance ratio. The LVG approximation is used for efficiently 
solving the radiative transfer equation, when the molecule level 
populations are not thermalized. LVG models are useful for a 
first order determination of the gas properties. However, as the 
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Fig. 7. Observed spectra towards the Overlap region (see Table 
IB. lb . See the caption of Fig. [T] 




NGC 4736 



n 



'^CO(3-2) 



"00(3-2) 
X 8 





200 400 

Velocity (km/s) 
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Markarian 231 



''C0(2-1) X 8 y 





C0(3-2) X 8 



Velocity (km/s) 



Fig. 9. Observed spectra towards Markarian 23 1 nucleus (see 
Table EB. See the caption of Fig. [T] 
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CENTAURUS A 
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Fig. 8. Observed spectra towards NGC 4736 nucleus (see Table 
IB.lt . See the caption of Fig. [T] 
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Fig. 10. Observed spectra towards Centaurus A nucleus (see 
Table EB. See the caption of Fig. [T] 



medium is assumed to be homogeneous, these models provide 
an average description of the molecular gas, which is known to 
exhibit structure at smaller spatial scales than sampled by these 
observations. 

Photo-dissociation region (PDR) models are more so- 
phisticated than LVG models as they solve simultaneously 
for the gas chemistry, photo-dissociation and thermal bal- 
ance, taking into account the relevant physical and chemical 
processes. Such models have been developed during the 
past two decades, for a variety of astrophysical sources, 
from giant molecular clouds illuminated by the interstellar 
radiation field to the conditions experienced by circum- 
stellar d isks or by dense mat ter , very clos e to hot massive 
stars jTielens & HollenbachI Il985allbl: iLe Bourlot et al 



19931: iKoester et al.l 11994 [Sternberg & Dalgarnol 11995 
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Fig. 11. Spectra of M 51 at various offset positions relative to 
the central position (0",0") given in Table [T] (see also Table 
IBTT i : (-24",-24") (top), (0",0") (middle) and (0",12") (bottom). 
See the caption of Fig. [1] The [CI]('^Pi-'^Po) sp e ctrum at the 



Fig. 12. Observed spectra towards M 83 nucleus (see Table 
IB. II) . See the caption of Fig. [T] 



central position (0",0") is from lOerin & PhilUpsI (l2000l) but it 
has been analyzed again to obtain an homogeneous dataset of 
observations. 



Hollenbach & TielensI 1 19991 iKaufman et al] Il999h . In FDR 



models, the line emission of the most abundant species, which 
determine the gas cooling rate, is usually calculated under the 
LVG approximation. Because of the many physical and chem- 
ical processes involved, PDR models use a simple geometry 
for the modelled cloud, which can be either plane parallel 
or spherical. When applied to CO line emission, PDR model 
predictions correspond to the integral along the line of sight 
of the contributions of the different regions in the modelled, 
plane parallel, cloud, with pronounced kinetic temperature and 
abundance variations from the warm outer shells to the cold 
interior Therefore, the contributions of these different regions 
to the CO emission depends on the rotational Une, the high-J 
CO lines being more sensitive to the warm gas (>20 K), while 
the CO(l-O) contribution is biased to the cold gas (<20 K). 

We have chosen to use both types of models as they are 
built with different hypotheses and therefore provide different 




Velocity (km/s) 



Fig. 13. Observed spectra towards Arp 220 nucleus (see Table 
IB. II ). Se e the caption of Fig. [1] In this figure, both spectra 
are from lOerin & PhilUpsI ( Il998h but they have been analyzed 
again to obtain an homogeneous dataset of observations. 
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Fig. 14. Observed spectra towards NGC 6946 nucleus (see 
Table IB. 11 1. See the caption of Fig. [T] 



information when the model outputs are compared with CO ob- 
servations. As LVG models are simple and fast, they are largely 
used in the analysis of extragalactic CO data. But LVG mod- 
els are obviously a crude approximation to the complex mix- 
ture of physical conditions in galaxy nuclei. Although solving 
the radiative transfer under the same LVG approximation, PDR 
models consider kinetic temperature and CO abundance gradi- 
ents along the line of sight. By comparing the PDR and LVG 
model predictions for the same source, we gain some insight 
on the quality of the models, especially on the reliability of the 
predictions for the missing CO lines, and on the computation 
of the CO cooling rate. 

5.2. LVG models 

5.2.1. Fitting procedures 

We have shown in lBavet et all (|2004 (Fig. 1) that the CO lines 
with J upper > 3 provide a good signature of the warm molecular 
gas {Tii > 20 K). In this paper, we used the same method, and 
constrain the fits of LVG or PDR models by the line intensity 
ratios from the observed CO submilhmeter lines. 



— 1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — 

, M83 arms 



C0(3-2) 



(90-100,30)- 



C0(3-2) X 10 




C0(4-3) 
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Fig. 15. Spectra taken in the M 83 arms at the offset positions 
(90"-100",30") (top) and (-90",-60") (bottom) (relative to the 
central position (0",0") listed in Table [1] see also Table IB. 11 1. 
See the caption of Fig. [T] 



More precisely, for IC 10 and IC 342, we used the follow- 
ing line intensity ratios (See Table ID : |,|:^^^~^| , ncod-l) ' 
TW^' "cQ(7-6) "hx)^-2) - observations of NGC 



4038 and Overlap suffer from large error bars. Indeed, the 
'^CO(4-3) and the '^CO(7-6) lines have a lower signal-to-noise 
ratio than other transitions (See Figs. [23] and l22l where these 
lines appear clearly too faint). Therefore, we chose for those 
sources the following ratios to constrain the models: t^^^-^. 



and 



'-CO(3-2) 



(See Table |4|. Since we missed the 



'2CO(4-3)' '2CO(6-5)' '2CO(6-5) 



C0(2-l) '-CO(3-2) 



"CO(3-2) 



'2CO(6-5) "CO(3-2} 

'2C0(7-6) line for M 83 and NGC 6946, we used the following 
line intensity ratios: "^'^-^^ "^'^-^^ 
(See Table ill. 

To compute these line intensity ratios, we used data identi- 
fied with an asterisk in Table |B.2| (preciselv. in the column 5 : 
the line area. A, in Kkms"') which have been previously scaled 
to a common beam size of 21.9"(see Sect. 14. 11 ). Values of line 
intensity ratios are listed in Table |4] The '^'CO data are partic- 
ularly useful for measuring the CO column densities while the 
'^CO data provide constraints on the kinetic temperature and 
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Table 3. Parameters of the "best fit" LVG and PDR models for IC 10, IC 342, NGC 4038, Overlap, M 83 and NGC 6946. For 
each galaxy, the same abundance ratio and the same Av" are used for constraining both LVG and PDR models. To compare 
observations and predictions deduced from models, the observed CO ('^CO and '^'CO) and the observed C cooling rates are listed 
at the bottom. 





IC 10 


IC 342 


NGC 4038 


Overlap 


M83 


NGC 6946 


13 r 


40 


40 


40 


40 


40 


40 


Av" (kms-') 


15 


54 


92 


155 


91 


158 


LVLr MUUbLb 














N( CO)/Ai' 

/ -2/1 -1\ 

(cm /kms ) 


3.0x10 


3.5x10 


3.2x10 


2.4x10 


6.0x10 


7.0x10 


n(H2)(cm-') 


7.0 X 10^ 


1.9 X 10' 


3.5x10' 


8.0x10' 


6.5 X 10' 


1.5x10' 




25 


40 


40 


145 


40 


130 




4.6 X 10""^ 


6.7 X 10"- 


1.7 X 10"- 


1.3 X 10"- 


3.6 X 10"-^ 


2.7 X 10"- 


LVG Predicted 
'-CO cooling 
(Wm"-sr"') 


5.7 X lO^' 


3.7 X 10"** 


4.1 X 10"** 


4.3 X 10"** 


1.3 X 10"' 


4.8 X 10"** 


LVG Predicted 
"CO cooling 
(Wm sr ) 


2.2 X lO^'" 


1.9 X 10"' 


1.0 X 10"' 


7.4 X 10"'° 


4.1 X 10"' 


1.0 X 10"' 


PDR MODELS 














FUV radiation 
(xGo) 


5.5 X 10* 


8.5 X 10' 


2.5 X 10= 


1.5 X 10' 


5.5 X 10'' 


1.5 X 10' 


n(H) (cm"') 


2.0 X 10' 


1.5 X lO*" 


3.0 X lO" 


3.5 X lO' 


5.0 X 10' 


1.0 X 10' 




1.3 X 10^^ 


2.9 X 10"- 


9.4 X 10"' 


9.2 X 10"' 


1.2 X 10"- 


2.9 X 10"- 


PDR Predicted 
'^CO cooling 
(Wm-^sr-') 


5.8 X 10"' 


3.6 X 10"** 


3.9 X 10"** 


4.9 X 10"** 


1.2 X 10"' 


4.6 X 10"** 


PDR Predicted 
'''CO cooling 
(Wm--sr-') 


3.1 X 10"'° 


3.2 X 10"' 


1.2 X 10"' 


1.8 X 10"' 


5.1 X 10"' 


2.1 X 10"' 


PDR Predicted 
C cooling 
(Wm--sr-') 


4.5 X 10"' 


2.8 X lO"** 


2.9 X lO"** 


2.6 X lO"** 


3.7 X 10"** 


8.7 X 10"** 


Observations'" 














Observed 
'^CO cooling 
(Wm-^sr-') 


4.0 X 10"' 


2.5 X lO"** 


1.7 X lO"** 


1.9 X lO"** 


4.3 X lO"** 


2.7 X 10"** 


Observed 
"CO cooling 

(Wm--sr-') 


5.7 X 10"" 


9.0 X 10"'" 


9.9 X 10"" 


1.2 X 10"'° 


4.2 X 10"'° 


3.9 X 10"'° 


Observed 
C cooling 
(Wm--sr-') 


1.1 X 10"' 


6.4 X 10"' 


9.2 X 10"'°'' 


7.1 X 10"'°'' 


2.3 X 10"''' 


1.4 X 10"''' 



" : deduced from Gaussian fits to the '^CO(3-2) line profiles. : Assumed tt^- ^ : FF is the filling factor of the CO emitting region in the 21.9" 
beam. '' : Based on the sole C(''Pi-^Po) transition. ; See C and CO lines identified with an asterisk in Table IB. 21 



the H2 density. The '^CO and '■'CO spectra used are displayed 
in Figs.[I]to [Tl(See Sect. lO ). 

We made use of a least square fitting method (taking into 
account the errors of the observed intensity line ratios) to deter- 
mine the physical conditions which reproduce the observations 



th e best. The line wid th Av and the abundance ratio (Xgak 



m 



axy 



Bavet et al.l (120041) ). are not part of the fitting process. Av is 



set to the value deduced from Gaussian fits of the '^CO(3-2) 



spectra (FWHM). The Av values used for the target galaxies 
are reported in Table |3] The abundance ratio is set to 40 
(see Table O for IC 10, IC 342, NGC 4038, Overlap, M 83 and 
NGC 6946, the value used in previo us LVG modelling works 
on the same galaxies (fo r IC 342 in iMauersbereer & Henkell 



1993l; lHenkel et al.lll998l: for NGC 4038 in lZhu et al.ll2003 



for 

M 83 inlMauersberger & Henkell ll993l and for NGC 6 946. in 
Israel & Baasll200lb . But for IC lO. IPetitpas & WilsonI (Il998al) 
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Table 4. Observed and predicted ("best" LVG model) Une ratios for the studied sources. 



A 


IC 10 


IC 342 


A 


NGC 4038 


Overlap 


A 


M83 


NGC 6946 


(Kkms-') 


obs.* 


obs.* 


(Kkms-') 


obs.* 


obs.* 


(Kkms-') 


obs.* 


obs.* 


'-CO(3-2) 


1.8±0.3'' 


2.2±0.3'' 


'-CO(3-2) 


1.5±0.2 


2.4+0.3 


'-CO(3-2) 


1.1+0.4" 


1.5+0.3" 


'-CO(4-3) 


'2CO(4-3) 


'-CO(4-3) 


"CO(3-2) 


2.5±0.5'' 


4.1 ±0.4° 


'-CO(3-2) 


1.1 ±0.2" 


1.7+0.8° 


'-CO(3-2) 


1.2+0.3" 


4.6+1.0" 


'^CO(6-5) 


'-CO(6-5) 


'-CO(6-5) 


'-CO(3-2) 
l-COf7-6) 


8.7±2.4" 


15.0±1.6° 


- 


- 


- 


- 


- 


- 


'-C0(2-l) 


10.8±4.6« 


18.5±5.3" 


'-C0(2-l) 


1.6+0.5° 


2.2+1.3° 


'-C0(2-l) 


1.4+0.1" 


6.9+1.6" 


'-CO(7-6) 


'-CO{6-5) 


l2CO(6-5) 


"CO(3-2) 


10.3±2.r 


5.5±0.7" 


'-CO(3-2) 


18.0+11.0° 


30.2+8.2°'' 


'-CO(3-2) 


8.5+2.6" 


17.1+2.0" 


^^^COf3-2) 






'-CO(l-O) 




11.3+0.4 


'-C0( 1-0) 






'-C0( 1-0) 




11.1+4.4 


'^COd-O) 


- 


"CO(l-O) 


- 


- 


'-'C0( 1-0) 


- 


'^CO(2-l) 


15.4±6.2 


10.9±2.6 


'-C0(2-l) 


24.3+9.7 


16.3+6.6 


'-C0(2-l) 


4.5+0.5 


29.7+6.7 


"C0(2-l) 


l-'C0(2-l) 


'-"^{2-1) 




LVG model 


LVG model 




LVG model 


LVG model 




LVG model 


LVG model 




In =25 K 


T^. =40 K 




=40 K 


=145 K 




T^, =40 K 


=130 K 


'^CO(3-2) 
12CO(4-3) 


1.2 


1.3 


'-CO(3-2) 
l-CO(4-3) 


1.1 


1.1 


'-CO(3-2) 
'-CO(4-3) 


1.1 


1.4 


'-CO(3-2) 


2.8 


4.1 


'-=CO(3-2) 


1.6 


2.7 


'-CO(3-2) 


1.4 


5.1 


12C0(6-5) 


12CO(6-5) 


l-CO(6-5) 


'^CO(3-2) 


8.4 


15.6 


'-CO(3-2) 


2.7 


5.9 


'-CO(3-2) 


1.9 


13.9 


^-C0(1-S) 


'2C0(7-6) 


'-CO(7-6) 


"C0(2-l) 


9.1 


17.9 


'-C0(2-l) 


1.5 


2.3 


'-C0(2-l) 


1.4 


6.1 




'-CO(6-5) 


'-CO{6-5) 


'-CO(3-2) 


11.6 


5.4 


'-CO(3-2) 


16.7 


23.9 


'-CO(3-2) 


10.1 


17.1 


13C0(3-2) 


l3CO(3-2) 


'3CO(3-2) 


'^CO(l-O) 
13CO(1-0) 


27.3 


4.2 


'-C0( 1-0) 
1 -'C0( 1-0) 


33.7 


36.8 


'-C0( 1-0) 
' 3 C0( 1-0) 


28.6 


18.6 


'■=^(2-1) 


14.4 


3.5 


'-C0(2-l) 


22.3 


27.9 


'-C0(2-l) 


14.8 


12.7 


13CO(2-l) 


l-'C0(2-l) 


'3CO(2-l) 



* : ratio derived from observations marked with asterisks in Table IB. 21 ° : values used as constraints for the LVG models; * : value from Table 
2 in lzhu et al., (,2003i) . 



used a abundance ratio of 50 and for Overlap, IZhu et al 



(|2003|) used two different abundance ratios (40 and 60). 



5.2.2. Results 

For these six sources, we varied N('^CO) from l.OxlO'* cm"^ 
to l.OxlO^o cm-2, Tk from 10 K to 255 K and n(H2) from 
10 cm-^ to 10^ cm . Model solutions (physical parameters 
which reproduce the observations the best) for each source 
are given in Table [3] As an example of the used fitting pro- 
cess, we present in Fig. [19] results of the LVG model cal- 
culations for NGC 4038 (right side) and M 83 (left side). In 
this figure, we plot Ax^ (x~'X^ ) contours in the 2D parame- 
ter space (n(H2)(cm"^) vs. Tic{K)). In this figure, the parame- 
ter N('^CO)/Av has been set to its "best fit" value. More pre- 
cisely, in Fig. [191 we have selected all LVG model solutions 
within a small interval in N('^CO)/Av centered around the best 
N('^CO)/Av value (see the chosen interval values at the top of 
each plot) for a better gridding of the (n(H2),Tk) parameter 
space. Similar plots were obtained for each target. The "best 
fit" model (with y- . ), is located at the intersection of the two 
black lines in Fig. [I9l Predicted line intensity ratios from the 
best fit models are listed in Table [4] In Fig. [T9]the best mod- 
els are located in the black zones which have a "banana" shape 
in the (n(H2), Tk) plane. It is clear that the fitting process us- 
ing the LVG framework is highly degenerated since acceptable 



fits (black zones) can be obtained over a large domain of the pa- 
rameter space. For NGC 4038, we conclude that the gas density 
is not well constrained while the temperature range is reason- 
ably narrow. For M 83, we have the reverse case: the gas den- 
sity is well constrained while the temperature range is broad. 
We notice that, despite the degeneracy, predicted CO line in- 
tensities computed from models located in the black areas in 
Figs. [19] do not differ significantly (variations in Une inten- 
sity between models localized in the black areas are < 3%). 
Therefore, LVG predictions are relevant for computing the CO 
cooHng rate. In Figs. [20] [IT] ISS] |23l |24]and |25l we present 
the best LVG models for the studied sources (IC 10, IC 342, 
NGC 4038, Overlap, M 83 and NGC 6946). For each source, 
we list the physical parameters corresponding to the best solu- 
tions in Table [3] Comparisons between observed and modelled 
line intensity ratios are given in Table [4] 

The LVG models have been used for predicting '^CO line 
intensities from '^00(1-0) up to '-CO(15-14) and the '^CO 
fine intensities from '^CO(l-O) up to '^CO(6-5). The predicted 
'^CO and the '^CO cooling rates are computed by summing all 
'^CO and '^CO line intensities and are listed in Table [3] We 
obtained the following results : 

* The "best" model for IC 10 is illustrated with grey triangles 
in Fig. [20] The lines which contribute the most to the '^CO 
and to the '^^CO cooUng rates are the ^^€0(5-4) line (26.6% 
of the total intensity), the '^00(6-5) line (25.1%) and the 
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Fig. 19. LVG model results for NGC 4038 (right side) and M 83 (left side). For each source, the "best" fit N( ^^CO)/Av (indicated 
at the top) is found by minimizing In fact, to avoid gridding problems, we use all models within a small interval centered 
around the best N('^CO)/Av value, listed in Table [3] We plot Ax^ contours in the Tk (K) vs. n(li2){cm^^) space. The contour 
levels for M 83 and for NGC 4038 are indicated in the plots. The LVG models with higher Ay- (>4) are not represented. In both 
figures, the dark areas represent the zones where is the lowest. Each cross corresponds to a LVG model solution. The best fit 
is located at the intersection of the two black lines. 



i2C0(4-3) Une (17.9%), the 1^00(4-3) Une (30.0%), the 
i3C0(5-4) line (29.5%) and the ^^€0(6-5) line (18.3%), 
respectively. 

* The "best" model for IC 342 is illustrated with grey 
triangles in Fig. |2T| The lines which contribute the most 
to the i^co cooling rate are '-CO(5-4) (28.1% of the 
total intensity), '2cO(6-5) (22.2%) and ^^€0(4-3) (21.7%) 
while the lines which dominate the '^CO cooling rate are 
^3CO(3-2) (36.3%), i^CO(4-3) (27.7%) and '^^CO(2-l) 
(19.2%). 

* For NGC 4038, we obtained a "best" model represented 
in Fig. |22]with grey triangles. The lines which contribute 
the most to the '^CO cooling rate are '-CO(6-5) line 
(22.8%), the '-CO(7-6) line (21.16%) and the i-CO(5-4) 
line (17.4%) while the main lines of the '^CO cooling rate 
are the i3CO(5-4) (33.0%), i^CO(6-5) (32.3%) and the 
'3CO(4-3) (22.9%) Unes. 

* For Overlap, we obtained a "best" model represented in 
Fig. |23]with grey triangles. The '^CO(6-5) line represents 
21.2% of the '^CO cooling rate while the '^00(5-4) and 
the ^2CO(7-6) lines correspond to 21.0% and 15.7%, 
respectively. The main lines for the '^CO cooling are the 



i-^CO(5-4) (29.9%), i^CO(4-3) (28.3%) and the i3CO(6-5) 
(22.9%) lines. 

* For M 83, the "best" model is represented in Fig. [24lwith 
grey triangles. The '^CO(7-6) line represents 22.7% of the 
i^CO cooling rate while the i"CO(6-5) and the '2cO(8-7) 
lines correspond to 19.4% and 18.4%, respectively. The 
'^CO cooling rate deduced from this "best" model is 
less accurate than for the previous galaxies since we 
constrained LVG models without the '^CO(7-6) line. The 
'^^CO cooling rate is dominated by the intensities of the 
i^^CO(6-5) (34.0%), i3CO(5-4) (32.8%) and the '3CO(4-3) 
(22.0%) lines. 

* For NGC 6946, we obtained a "best" model represented in 
Fig. [25] with grey triangles. The main lines for the '-CO 
cooUng rate are the i2cO(5-4) (23.6%), the i2cO(4-3) 
(20.3%) and the i-CO(6-5) (18.9%) lines. The '^CO cool- 
ing rate deduced from this "best" model is less accurate 
than for the previous galaxies since we constrained LVG 
models without the '^CO(7-6) line. The '^CO cooling rate 
is dominated by the intensities of the '''CO(3-2) (28.8%), 
i^CO(4-3) (27.5%) and the i^CO(5-4) (18.4%) lines. 
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Fig. 16. Spectra taken towards the NGC 6946 arms 
positions (llO'MOO") (top) and (150",-20") (botto 
to the central position (0",0") listed in Table [T] ses 
IB. lb . See the caption of Fig. [1] 



We present in Figs. |20l EU |22l |23] EH and |25] 
the predicted integrated areas (A in Kkms ') and the pre- 
dicted Une intensities (1 in Wm"^ sr ') for all '^CO transi- 
tions up to J=15-14 and for all '''CO transitions up to '''CO(6- 
5). We have computed the predicted integrated areas. A, us- 
ing the predicted antenna temperature from the models and 
multiplying it by the line width Av and by the surface fill- 
ing factor, FFlvg (both listed in Table O. For each ob- 
served CO line, we have estimated the surface filling factor 
(FFlvg) of molecular clouds in the 21.9" beam using the ratio 
Aofcjpn'fl//o«s(Kkms"')/A„,o^/e/(Kkms"'). We have made an aver- 
age, weighted by the S/N, for the filling factors pertaining to 
each observed transition. For NGC 4038 and Overlap, we have 
not taken into account the '^CO(4-3) and the '-CO(7-6) lines 
because they have the lowest signal-to-noise ratio (see Figs. |22] 
and l23T l. For the studied sources, we have obtained surface fill- 
ing factor of a few % (1.3 to 6.7%) in a 21.9" beam size. 




10 -10 
RA offset (arcsec) 



Fig. 17. Map of the velocity integrated intensity of '-CO(3-2) 
(top) and '2cO(6-5) (bottom) towards IC 342. The CO emis- 
sion is integrated over the velocity ranges - 200 to 200 kms"' 
and -50 to 150 kms"', respectively. For both maps, the crosses 
show the observed positions. The contours of the '^CO(3-2) 
map range from J T^bdv - 12 Kkms"' to 1 12 Kkms"' with 10 
Kkms"' steps. For the '^CO(6-5) map, the contours range from 
jT,nbdv = 8 Kkms"' to 108 Kkms"' with 10 Kkms"' steps. 
The intensity peak value of the '^CO(3-2) map is J T„,hdv ~ 
110 Kkms"'. The intensity peak value of the '^CO(6-5) map 
is j T,„t,dv X 108 Kkms"'. The black box in the upper figure 
represents the size of the lower figure. The circles indicate the 
beam sizes at the frequency of the '^CO(3-2) line (top) and at 
the frequency of the '^CO(6-5) line (bottom). No pointing error 
was detected for this map. 



5.2.3. Discussion 

From the above analysis we can conclude that LVG models can 
be confidently used to determine the CO cooUng rates, pro- 
vided that a sufficiently complete data set is available. Even 
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ing rates. We also showed that '^CO lines with J„pper > 
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Fig. 18. Map of the velocity integrated intensity of '^CO(3-2) 
(top) and '-CO(6-5) (bottom) lines towards M 83. For both 
maps, the CO emission is integrated over the velocity range 
200 to 800 kms \ and the crosses show the observed positions. 
The contours of the '^CO(3-2) map range from J T„,hdv - 3 
Kkms"' to 18 Kkms"' with 3 Kkms"' steps, and from J T„,i,dv 
= 21 Kkms"' to 102 Kkms"' with 10 Kkms"' steps. For the 
'^CO(6-5) map, the contours range from J T,„i,dv - 20 Kkms"' 
to 320 Kkms ' with 25 Kkms"' steps. The intensity peak value 
of the '2cO(3-2) map is J T„,hdv ^ 100 Kkms^'. The intensity 
peak value of the '^CO(6-5) map is j T,„i,dv a; 318 Kkms"'. 
The '~CO(6-5) map has been shifted by +5" along the minor 
axis due to a pointing error The black box in the upper figure 
represents the size of the lower figure. The circles indicate the 
beam size at the frequency of the '^CO(3-2) line (top) and at 
the frequency of the '^CO(6-5) line (bottom). 



though several points in the parameter space with different 
physical conditions fit the observed data equally well, the de- 
rived CO cooling rates are very similar in all studied cases. 
The key point is the combination of '^CO(6-5) and '^CO(7- 
6) data as these two lines provide the largest contribution to 
the observed '^CO cooling. It will be particularly interesting 
to observe the '^CO(5-4) line which is shown to be also a 
main cooling line in these galaxies. In addition to the '■'CO(5- 
4), i2cO(6-5) and i-CO(7-6) lines, data for '"CO(8-7) and 
'^CO(9-8) lines would be most useful to discriminate mod- 
els, and for a more accurate determination of the '^CO cool- 



are predicted to be weak and will not have significant antenna 
temperatures (see Figs. |20l E] |22] |23l Ell and |25] plots 
on the left side). In addition '^CO(6-5) data would also be ex- 
tremely useful for constraining the models, for better estimat- 
ing the '^CO cooling rates and for measuring the opacity of 
'^CO(6-5) line. Indeed, when comparing the predicted and ob- 
served intensities of the '^CO lines, the importance of high-J 
'^CO lines shows up. For the studied galaxies, ' ^CO(3-2) is 



the m ost intense observed line (see this work and lBavet et al 
|2004 . For most sources, predicted '-''CO(6-5) line intensities 
are at least as strong as '■'CO(3-2) intensities. 

Because the LVG models have been constrained by the 
high-J CO line intensity ratios, intensities of the low-J CO tran- 
sitions are not well fitted (see for instance the '^CO(l-O) or the 
'^CO(2-l) predicted intensities in Figs.l20l 1211 [22]and l25T l. 
To relieve that problem, a two component LVG model would 
be needed, one component fitting the low-J CO transitions and 
another one f or the high-J CO tran s itions. It was not our pur- 
pose here but lHarrison et al.l (1 1999b : [Bradford et aD (l2003b did 
it for the galaxy NGC 253. 

We have compared the physical parameters (See Table 
|3]l corresponding to our best fit LVG models with results 
obt ained in previous s tudies for the same galaxies. For IC 
10, IPetitpas & Wilson ( 1998a ) presented an acceptable LVG 



solution: N('-'CO)/Av =5.0x10" cm-^/kms"', n(H2)= 10"*- 
10^ cm"^ and Tk -lOQ K. We suggest a higher value for 
the gas density and lower values for N('^CO)/Av = 3 x 
10'^cm~^/kms~' and Tk - 25 K than those proposed by 



Petitpas & WilsonI (Il998al) . The differences can be explained 



by the fact t hat we do not use t he sam e set of CO line inten- 
sity ratios : iPetitpas & WilsonI (Il998al) used line ratios com- 
bining low-J CO lines (12C0(2-1), '^00(2-1), '^00(3-2) and 
'^CO(3-2)) while we used a larger number of intensity ra- 
tios focussed on th e high -J CO lines (See Table |4]i. For IC 



342, llsrael & BaasI (120031) fit their data using a LVG model 
with N(i2cO)/Av =6-10x10 '^ cm'^/kms ^ ', T;^ =100-150 K 
and n(H2)= 3.0x10^ 
model with N('2cO)=3-4xlO 
sity around njUj) - 2.0x10^ 



Eckart et al. d 19901) suggested a 



cm , Tk >20 K and a den- 
Q"-^ for the center of IC 342. 



Meier et al. deduced from their observations a beam- 

averaged density of n(H2)= 1.3x10'' cm Our model agrees 
very well with the gas density values proposed in the two latter 
articles. T he fitted N('^CO)/Av valu e is intermed i ate bet ween 



results of llsrael & Baal (l2003l) and lEckart et al.1 (Il990l). For 



Tif, w e agree better with the value proposed by Eckart e t al 
d 19901) than with the one prop osed by llsrael & Baas (2003). For 
NGC 4038. IZhu et al.l ( l2003b proposed a one component model 



(N('2C0)/Av =3.4x10'^ cm-2/kms-i, n(H2)= 1-4x10^ cm" 
and Tk =43 K) in good agreement with our study except for the 
n( H2) value. For Ov erlap, the one component model proposed 



bv lZhuetal.1 1003 has NC^CO)/Av ^1.0x10'" cm^-'/kms^ 



n(H2)= 1-5x10^ cm-^ and Tk =33 K for their (0":10") ofi^- 
set positions which almost corresponds to our Overlap po- 
sition. The agreement of the former study with our values 



14 



is not as good as for NGC 4038 but IZhu et al.l (l2003k used 
a higher CO abundance ratio ('^CO/'^CO = 60) than we 
do (4 0), which may be at the origin of the differences. For 
M 83, llsrael & BaasI ( 1200 lb presented two LVG models: one 
with N(i2cO)/Av =1-3x10'^ cm~^/kms-\ Tk =30-150 K and 
n(H2)=0.5-3.0xlO^ cwT^, and one with ^C^CO)IAv =0.06-1 
xlO'^ cm--/kms-i, Tk =60-100 K and n(H2)= 0.03-1.0x10^ 
cm"^. Our fitted parameters compare well with both mod- 
els, although the y are closer t o the second model conditions. 
For NGC 6946. iwalsh et all (l2002h proposed the following 
LVG model: N(i2cO)/Av =2.9x10'^ cm;^/km s-', Tk =40 K 



and n(H2)= 2.0x10^^ cm~^. llsrael & BaasI (120011) presented two 
other possible LVG models: one with N('^CO)/Av = 1 - 10 x 
10'^ cm--/kms-', = 30 - 150 K and n(H2)= 0.5 - 1.0 x 10^ 
cm-^ and one with N('-CO)/Av = 3 - 6 x 10^'' cm~2/kms"\ 
= 30 - 150 K and n(H2)= 0.1 - 1.0 x 10"* cm'^. All 
these models are consistent with the parameters obtained in the 
present work. 



5.3. PDR models 
5.3.1. Fitting procedures 

We use the PDR models deve loped by Le Bourlot et al. (1993) 
for Galactic sources (see also lLe Petit et al.ll2002l) . The source 
is modelled as a plane-parallel slab, illuminated on both sides 
by FUV radiation to better reproduce a starburst environment 
where massive stars, and giant molecular clouds, are spatially 
correlated. Model parameters include the gas density, assumed 
uniform, the intensity of the illuminating FUV radiation, the 
gas phase elemental abundances, the grain properties and the 
gas to dust ratio. 

We considered that all galaxies have a metallicity close to 
solar (See Table [TJ, and we used Milk y Way abundances in 
our models (12 + log(§) =8.90 ± 0.04 lArimoto et al.lll996h . 
For IC 10 and IC 342, this hypothesis may be rather crude con- 
sidering the lower and higher (respect ively) metallicity of t heir 
nuclei (12 + /og(g) = 8.31 ± 0.2 inkrimoto et al.lll996l and 
l2 + logi§) =9.30 in lWa-Costas & Edmundsll992h Therefore 
results obtained on these sources should be considered more 
cautiously. Standard grain properties and a gas to dust ratio 
appropriate for Galactic interstellar clouds have been adopted. 
The '^C/'^'C ratios are the same as used above for the LVG 
models (see Table O. We have sampled a wide range of pa- 
rameter space, varying the gas density, n(H) from 1.0x10^ 



cm--^ to 1.0x10^ cm , and the incident FUV flux, xfuv, from 
1.5 X lO-' Go to 1.5 X 10'' Go, where Go is the local average 
interstellar radiation field (ISRF) determined by Draine (1978) 
(Go = 2.7 X 10"^ ergcm--s-^). We fit the PDR models with 
the same line intensity ratios used to fit the the LVG models 
(See Table [5]l but computed from values listed in Table IB. 21 
(with an asterisk) expressed in Wm"-sr"', a unit more appro- 
priate for PDR model calculations. Indeed, the output values 
(emissivities) from PDR models calculations are expressed in 
ergcm^-s^'sr"'. 



The best model is obtained using a least square fitting 
method as applied in Sect. 15.21 taking into account the errors 
on the observations. 

5.3.2. Results 

Best fit solutions (physical parameters) for the selected galax- 
ies are presented in Table [3] Model predictions for the CO line 
emissivities are shown in the right side of Figs. [20l |2T1 |22l 
|23l |24]and |25]for IC 10, IC 342, NGC 4038, Overlap, M 
83 and NGC 6946, respectively. Model predictions have been 
scaled by a surface filling factor (FFp^R listed in Table [3]) to 
match the observed line intensities. As stated above, PDR mod- 
els have been developed for local interstellar clouds. The veloc- 
ity dispersion of the modelled cloud is a parameter in those 
models, which is used for computing the photo-dissociation 
rates of H2 and CO, and the line emissivities. This parame- 
ter i s set to 1 kms~\ a ty pical figure for local molecular clouds 
(see lWolfireetalJll990l) . Using a PDR model for fitting galaxy 
observations is complicated by the fact that, in a galaxy, many 
PDRs contribute to the signal detected in each beam, result- 
ing in a broad line (tens to hundreds of kms"'), compared to a 
single PDR line (1 kms"'). To correct for this effect, the mod- 
elled emissivities have been multiplied by the line width ratio, 
^hv(PDR) ■> where lS.v{galaxy) is the line width of the CO lines 
as reported in Table [H and where Av{PDR) =1 kms Once 
this correction is performed, PDR model results are compared 
with observed data in the same way as the LVG models (see 
Sect. 15.21 ). The surface filling factors of the emission in the 
beam, FFpdr, are also computed (and listed in Table O and 
they are in reasonable agreement with those derived from LVG 
models for all sources. 

The PDR solutions obtained from the least square fitting 
procedure are described below for the individual galaxies. The 
'^CO and the '^^CO cooling rates are summarized in Table [3] 

* For IC 10, the "best" model is represented in Fig. [20] with 
grey squares. The '^CO(5-4) line represents 24.7% of the 
total '-CO cooHng rate while the '^00(6-5) and '^00(4-3) 
lines correspond to 24.6% and 16.8%, respectively. The 
main lines for the '^CO cooling rate are the '■^CO(5-4) 
(33.0%), 1^00(6-5) (26.8%) and '^00(4-3) (25.8%) lines. 

* For IC 342, the "best" model is represented in Fig. |2T| 
with grey squares. The lines which contribute the most 
to flie '^CO cooling rate are the '-CO(5-4) (27.3% of the 
total intensity), '^00(6-5) (23.6%) and '^00(4-3) (20.1%) 
lines. The lines which contribute the most to the '^CO 
cooling rate are the '3CO(4-3) (30.6%), '^CO(5-4) (29.6%) 
and 1^00(3-2) (17.7%) lines. 

* For NGC 4038, the "best" model is represented in Fig. |22] 
with grey squares. The lines which contribute the most to 
the '-CO and to the '^^CO cooling rates are the '2CO(6-5) 
(23.4% of the total intensity), 1^00(5-4) (20.3%), ^^COO- 
6) (19.0%) lines and the '^^CO(6-5) (35.7%), 1^00(5-4) 
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J(upper) J(upper) 

Fig. 20. "best" LVG and PDR models compared with observations for IC 10 (See Table [3] for the model parameters) . On the 
left side, we plot integrated intensities in Kkms ' vs iupper- On the right side, we plot I in Wm^^sr ' vs iupper- In all figures, grey 
triangles represent the LVG model while the grey squares represent results from the PDR model (see Sect. I5.3l l. Observations 
(with error bars) taken from the literature and from our data set are shown in black. The '^'CO transitions are represented with 
open symbols while the '^CO transitions use filled symbols. 
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Fig. 21. "besf'LVG and PDR models compared with observations for IC 342 (See Table [3]for the model parameters). See caption 

r>f Birr 
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Fig. 23. "best" LVG and PDR models compared with observations for Overlap (See Table |3]for the model parameters). See 
caption of Fig. |20] 
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Fig. 24. "best" LVG and PDR models compared with observations for M 83 (See Table [3]for the model parameters). See caption 
ofFig.EOl 
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Fig. 25. "best" LVG and PDR models compared with observations for NGC 6946 (See Table |3]for the model parameters). See 
caption of Fig. |20] 



18 



(33.4%) and '^€0(4-3) (21.1%) lines, respectively. 

* For the Overlap region, the "best" model is represented 
in Fig. |23]with grey squares. The lines which contribute 
the most to the '^CO and to the '^CO cooling rates are 
the '2CO(6-5) (22.9% of the total intensity), ^^CO(l-6) 
(19.7%), '2CO(5-4) (19.1%) lines and the '^^€0(6-5) 
(34.7%), '^CO(5-4) (33.3%) and '^€0(4-3) (21.6%) Unes, 
respectively. 

* For M 83, the "best" model is represented in Fig. |24] 
with grey squares. The lines which contribute the most 
to the '^CO cooling rate are the '^00(6-5) (20.4% of the 
total intensity), ^^CO{7-6) (20.4%) and i2cO(5-4) (15.4%) 
lines. The lines which contribute the most to the '^CO 
cooling rate are the '^€0(6-5) (34.1%), '^€0(5-4) (32.9%) 
and i3CO(4-3) (21.9%) lines. 

* For NGC 6946, the "best" model is represented in Fig. |25] 
with grey squares. The lines which contribute the most 
to the '^CO cooling rate are the '2cO(5-4) (28.5% of the 
total intensity), '2cO(4-3) (22.7%) and '2cO(6-5) (21.6%) 
lines. The main lines for the '^CO cooling rate are the 
^3CO(5-4) (33.0%), '3CO(4-3) (27.8%) and '^^CO(6-5) 
(23.7%) hnes. 

For the IC 10 nucleus (Fig.l20t. we obtained a good agree- 
ment between the two solutions proposed by the LVG and PDR 
models (similarly shaped curves and '^CO cooling rates); prob- 
ably because of the high data quality. Small differences appear 
in the '^CO cooling curves since the peak seems to be located 
at the '^CO(4-3) transition for the LVG model while the PDR 
peak clearly corresponds to the '^CO(5-4) transition. '^CO(6- 
5) observations would be very useful to discriminate between 
the model solutions. 

In IC 342 (Fig. |2TI ). we observed an excellent agreement 
between the two model solutions, the sole difference being in 
the peak position for the '^CO cooling curves (the LVG peak 
corresponding to the '^^CO(3-2) line while the PDR one being 
located at '^^CO(4-3)). Here also, '^CO(6-5) data would be very 
useful to better localize the '^CO peak. In this case, we have a 
very complete and high quality dataset. The '^CO cooling rate 
deduced from the PDR model is very close to the '^CO cooling 
rate computed from the LVG model. 

At the Antennae positions (Overlap and NGC 4038 in Figs. 
I22I and I23] ). LVG models, as well as PDR models, do not 
reproduce observations adequately due to the low signal-to- 
noise ratio of both '^00(7-6) and '-CO(4-3) Unes. Indeed, be- 
cause these latter observations are weaker than we expected, we 
have not considered them in the least square fitting procedures. 
However, for NGC 4038, both model solutions (LVG and PDR) 
are coherent : the maximum of the '~CO cooling curve deduced 
from the LVG model is the same as the result from the PDR 
model ('^CO(6-5)) and the '^CO cooling rates are similar. For 
Overlap, the maxima of the LVG and PDR '^CO cooling curves 



appear both for '^CO(6-5) and the predicted '^CO cooling rates 
from the LVG and PDR models are close to each other Despite 
this, we clearly see differences in the predicted line intensi- 
ties from LVG and PDR models, certainly due to the lack of 
relevant '-CO(7-6) detections. Concerning the '^CO cooling 
curves shown in the right side of Fig. |22](for NGC 4038) the 
two model solutions are compatible. Consequently, the total 
'^CO cooling rate agrees between the LVG and PDR models 
(see Table [3]). For the Overlap position (Fig. |23] ), the LVG 
and PDR '^CO cooling curves peak at the different positions 
(i^CO(6-5) for the PDR and '-^00(5-4) for the LVG) but the 
'^CO cooling rate are different by a factor < 3 (see Table O. 
For NGC 4038 and Overlap, '^^CO(6-5) observations would be 
very useful to discriminate between model solutions. 

For M 83, the lack of '-CO(7-6) data explains the differ- 
ences (noticeable in Figs. [24l i between the LVG and PDR so- 
lutions for the '^CO cooling curves. The peak corresponds to 
the '^00(7-6) line in the LVG and PDR models. Observations 
of the '-CO(7-6), '^00(8-7) Unes and possibly up would be 
very useful. Despite the different shape of the '^CO cooling 
curves, very similar '"'CO cooling curves are obtained. We no- 
tice that although the M 83 dataset is less complete than for 
other sources, we obtained fairly similar values of the LVG and 
PDR '^CO cooUng rates (see Table O. 

For NGC 6946, the peak positions of the '"CO cooling 
curves are similar ('2cO(5-4)) for the LVG and the PDR 
models. Also in this case, observations of the '^CO(7-6) and 
'~CO(8-7) lines appear essential to discriminate between model 
solutions. For the '-^CO cooling curve, we observed stronger 
differences especially for the peak position since the LVG max- 
imum corresponds to the '^CO(3-2) line while the PDR maxi- 
mum appears for the '^CO(5-4) line. Once again, the '^CO(6- 
5) observation would be very useful to better localize the '^'CO 
peak. 

5.3.3. Discussion 

As PDR models provide predictions for lines of C^, O, and 
C together with CO rotational lines, it is interesting to check 
whether the proposed PDR solutions are appropriate to repro- 
duce observation of other cooling lines such as the fine struc- 
ture lines of atomic carbon. In fact, for all sources studied here 
(see Table O, the predicted C cooling rates are typically larger 
than those observed one by significant factors (between 4.5 and 
50). 

For a large majority of the sample galaxies, we conclude 
that PDR models focussed on the high-J CO lines do not cor- 
rectly reproduce the observed atomic carbon data well. To bet- 
ter model the observed atomic carbon transitions, another set 
of physical conditions might be a better choice. For example, 
we could have increased the total extinction A,,, set to 10 mag 
in the PDR models we used. It is expected that the relative con- 
tributions of C and CO will change as CO lines are produced at 
larger depths than C. We performed a few calculations, variyng 
the total extinction across the slab. A,,. As expected, Co lines 
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Table 5. Observed and predicted ("best" PDR model) line intensity ratios. 



1 


IC 10 


IC 342 


1 


NGC 4038 


OVERLAP 


¥V 111 i3l J 


obs.* 


obs.* 


y YV 111 31 J 


obs.* 


obs.* 


'-CO(3-2) 


0.8±0.1'' 


0.9±0.1" 


'-CO(3-2) 


0.6±7.7xl0-2 


1.0±0.1 


■^00(4-3) 


'2CO(4-3) 


'-CO(3-2) 


0.3±6.3xl0-2'' 


0.5±4.5xl0-2° 


'-CO(3-2) 


0.1 ±2.4x10-'° 


0.2±9.4xl0--° 


'-CO(6-5) 


'-CO(6-5) 


'^CO(3-2) 


0.7±0.2'' 


1.2±0.1'' 


'-CO(3-2) 


1.8±0.6'' 


1.2±0.2° 


12CO(7-6) 


'2C0(7-6) 


'^C0(2-i) 


o.3±o.r 


0.4±0.1 


'-C0(2-l) 


(5.8±1.9)xlO-2'' 


(8.3±4.9)xlO-2-'' 


COO-6} 


CO(6-5) 


'-CO(3-2) 


11.8±2.3" 


6.3±0.8" 


'-CO(3-2) 


20.6±12.6'' 


34.5+9.4°' 


'3cO(3-2) 


l-'CO(3-2) 


'^CO(l-O) 




12.9±0.5 


'-C0( 1-0) 






"CO(l-O) 


- 


' 3 C0( 1-0) 


- 


- 


'-C0(2-l) 


17.6±7.0 


12.5±3.0 


'-C0(2-l) 


27.9+11.1 


18.6+7.5 


'3CO(2-l) 


'3CO(2-l) 




PDR model 


PDR model 




PDR model 


PDR model 


'-CO(3-2) 


0.5 


0.5 


'-CO(3-2) 


0.4 


0.4 


12CO(4-3) 


'-CO(4-3) 


'^CO(3-2) 


0.3 


0.4 


'-=CO(3-2) 


0.2 


0.2 


12CO(6-5) 


12CO(6-5) 


'-CO(3-2) 


0.5 


0.9 


'-CO(3-2) 


0.3 


0.2 


l-CO(7-6) 


l-CO(7-6) 


'-C0(2-l) 


0.1 


0.3 


'-C0(2-l) 


4.5x10-- 


4.9x10-- 


l2CO(7-6) 


'-CO(6-5} 


i-CO(3-2) 


12.3 


6.3 


i-CO(3-2) 


19.6 


15.7 


13C0(3-2) 


'3CO(3-2) 


'-C0( 1-0) 


18.9 


11.3 


'-C0( 1-0) 


24.6 


22.6 


' 3 C0( 1-0) 


1 -'C0( 1-0) 


'^C0(2-l) 


14.6 


7.1 


'-C0(2-l) 


21.9 


18.5 


'3CO(2-l) 


'3C0(2-1) 



■'CO(3-2) 



■CO(4-3) 



-CO(3-2) 



-CO(6-5) 
■^CO(3-2) 



■CO(6-5) 



CQ(3-2) 



'C0(3- 



C0( 1-0) 



3CQ(l-0) 



'C0(2-l) 



-<C0(2-1) 



C0{3- 



-CO(4-3) 



C0(3- 



-CO(6-5) 
■^CO(3-2) 



■con -6) 



C0(2-l) 



-CQ(6-5) 



■C0(3- 



3CO(3-2) 



CO( 1-0) 



3CO(1-0) 



■'C0(2-l) 



'CO(2-l) 



M83 
obs.* 



0.5+0.2° 



0.1+3.9x10- 



(5.2+3.8)xlO- 



32.7+10.1° 



5.1+0.6x10-- 



PDR model 



0.4 



0.2 



0.2 



5.1x10-- 



9.8 



19.3 



12.6 



NGC 6946 
obs.* 



0.6+0.1° 



0.6+0.1° 



0.3+5.9x10- 



19.5+2.3° 



12.7+5.1 



34.0+7.6 



PDR model 



0.5 



0.5 



1.0 



0.1 



18.3 



22.2 



20.1 



* : r atio derived from observations marked with asterisks in Table IB. 21 ° : values used as constraints for the LVG models; * : value from Table 
2 in lZhu et al ; values not observed neither found in literature. 



become more prominent relative to C lines with increased A,.. 
However, the observed ratio is never reached. A more com- 
pleted study of the parameter space is needed, which is beyond 
the scope of this paper. 



6. Comparison with the Center of IVIilky Way and 
the Cloverleaf QSO 

It is interesting to compare this galaxy sample with other, well 
known, sources. As we previously discussed previously, the 
peak position of the '"CO cooling curves shifts between nuclei. 
In Fig. |26]and |27] we compare the '^CO cooling curves for 
the center of the Milky Way and the Cloverleaf QSO [HB89] 
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1413+117 with our sample galaxies. We completed the set of 
observed '"CO lines with predictions from PDR models up to 
'^CO(15-14). Because the PDR model for Overlap is less re- 
liable, we excluded this source from this study. We have kept 
NGC 4038 because LVG and PDR model results are consis- 
tent for this object. Flux values for the center of the Milky Way 
and for the Cloverlea f QSO are from Fixsen et al. (19§9) an d 



1999h : IWeiB et al.l (120031) . 



Barvainis et al.l(ll997l) : lTsuboi et al 
respectively. 

The CO cooling rates (listed in Table O appear to be higher 
in the starburst environments (NGC 253 in iBavet et al. 2004 , 
M 83) than in the normal spiral gal axy (IC 342 NGC 6946) 
or irregular galaxies (Henize 2-10 in lBavet et alj|2004l and IC 
10). This phenomenon seems not to be solely a distance effect 
since we obtained the highest observed CO cooling rates for 
NGC 253 and M 83, two galaxies which are neither the nearest 
nor the farthest sources (2.5 Mpc and 3.5 Mpc, respectively). 
Moreover, the lowest observed CO cooling rates are obtained 
for Henize 2-10 and IC 10 which distances are 6 Mpc and 1 
Mpc, respectively. 

By looking at Fig. |26]and |27] we can distinguish three dif- 
ferent behaviors ; for Henize 2-10 and the Cloverleaf QSO, the 
CO cooling curve peaks at the '^CO(7-6) line or up; for IC 10, 
NGC 253, NGC 4038, M 83 and NGC 6946, the turnover ap- 
pears at the '^CO(6-5) lin^ and for the center of the Milky 
Way and IC 342, the turnover is found near the '^CO(4-3) 
line. We may ask whether these observed differences are solely 
consequences of the difference in linear resolution solely, or 
whether they are due to differences in physical conditions. 

The turnover of the '-CO cooling curve does not depend 
on the distance, since there is no obvious correlation between 
the distance and the position of the CO peak. For instance, the 
fact that the Cloverleaf QSO and Henize 2-10 have very similar 
CO cooling curve shapes (See Fig. l26b , is related to the close 
similarities of the physical properties of the warm gas for these 
two sources, which translates into similar CO line ratios. 

The position of the peak in the CO cooling curves is a par- 
ticularly interesting parameter, as it moves towards higher J for 
more actively star forming galaxies. Indeed, such curves could 
be used to diagnose the gas properties and measure the cooling 
rate in galaxies. Few detections are necessary (high-J CO lines 
as '^CO(6-5) and '^CO(7-6)) to obtain a first estimation of the 
gas cooling rate. 

7. Conclusions 



In this paper, we have presented observations of the C and CO 
submillimeter lines (up to the '^CO(7-6) line at 806 GHz and 
the [CI](^P2-^Pi) Une at 809 GHz) for the Antennae galaxies 
(NGC 4038 and the overlap region between the two nuclei). 



- For NGC 6946, due to the lack of the '-CO(7-6) detection, the 
LVG and PDR predictions are less reliable than for the other sources 
for which we have a more complete dataset. To determine the position 
of the NGC 6946 CO cooling curve, we chose to trust the observations 
rather than the predictions. 
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Fig. 26. Flux (Wm ^) vs. Jupper for the center of the Milky 
Way (top left, triangles), IC 342 (top left, points), NGC 6946 
(top right, points), M 83 (top right, triangles), the Cloverleaf 
QSO (bottom left, triangles), Henize 2-10 (bottom left, points) 
IC 10 (bottom right, triangles) and NGC 253 (bottom right, 
points). PDR models (grey points or grey triangles) have been 
used to obtained line fluxes of unobserved lines (eg. '^CO(5-4), 
'2CO(7-6), '-CO(8-7), "200(9-8),...) depending on the com- 
pleteness of the dataset obtained for each source. Observations 
are shown in black with error bars. To ease the comparison, 
we applied scaling factors on all sources data except for NGC 
253. The value of these scaling factors is specified for each 
source on the plots. We notice that for th e Cloverleaf QSO, t he 
'2CO(1-0) line flux is an upper limit fsee lTsuboi et al.lll999h . 



the nuclei of IC 10, IC 342, M 83 and NGC 6946 (see Tables 
IB. H and |B.2| l. We also detected submillimeter C and CO lines 
for the following galaxies: nucleus of Arp 220, Centaurus A, 
IRAS 10565+2448, M 51, M 82, Markarian 231, NGC 3079, 
NGC 4736, NGC 6090 nuclei and in the spiral arms of M 83 
and NGC 6946. 
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Fig. 27. Flux (Wm--) vs. 7„^,,,, for the NGC 4038 nucleus 
(points). See the caption of Fig.[26l In this figure, we plot ob- 
servations of the '^CO(4-3) and '-CO(7-6) lines even though 
they have a low signal-to-noise ratio. We also plot their PDR 
predicted fluxes. 



We succeed in fitting all observed CO lines to accurately 
estimate the CO cooling rate. LVG and PDR models have been 
used for fitting the CO line intensities. As we selected sources 
with different morphological types, we have compared the con- 
tribution of C and CO and the shape of the CO cooling curves 
between galaxies. 

1 The total CO cooling rates were estimated summing 
intensities (I in Wm "sr ') of all CO transitions up 
to '-CO(15-14), derived either from observations or 
from PDR model predictions. We showed that the total 
CO cooling rate is higher in the starburst nuclei (M 
83 and NGC 253) than in normal galaxies (IC 342 and 
NGC 6946) and irregular galaxies (IC 10 and Henize 2-10). 



83 and NGC 6946. As we missed the [CI](^P2-^^Pi) line at 
809 GHz for the Antennae, M 83 and NGC 6946, the latter 
figures should be confirmed with more observations. 

Our analysis made use of two complementary models: LVG 
models (Large Velocity Gradient for radiative transfer calcula- 
tions) and PDR models (combining chemistry, thermal balance 
and radiative transfer) for deriving the physical conditions in 
the galaxy nuclei, and for computing the intensity of the miss- 
ing CO lines. We obtained the following results: 

4 We succeeded in deriving physical properties (the kinetic 
temperature (T^:), the gas density (n(H2)), the CO column 
density divided by the line width (N('^CO)/Av) and the 
Far-UV radiation field (xfuv)) of the warm and dense 
molecular gas, using high-J CO and '^CO line intensity 
ratios for constraining the models. The LVG models alone 
can be highly degenerated as an extended region in the 
(n(H2); Tk) parameter space provides equally good fits. 
However, best models provide consistent estimations of 
the CO cooling rates. 

5 The surface filling factor of the molecular emission 
amounts to 1-7% in a 21.9"beam. High resolution obser- 
vations, using existing (IRAM-PdBI, SMA) or planned 
(CARMA, ALMA) (sub)minimeter interferometers, 
should be able to resolve individual sources in these 
galaxies. 

6 The '^CO(5-4) Une appears as one of the main CO cooling 
line for most sources. Unfortunately, the line is blocked 
by the atmosphere for local galaxies. Spatial missions 
(eg the Herschel Space Observatory) will soon provide 
measurements of this important line. We also showed 
the key role of the '3CO(6-5) and ^-CO{8-7) Hnes as 
constraints to PDR and LVG models. Future observations 
with ALMA and the Herschel satellite will provide these 
informations. 



2 The shape of the CO cooling curve depends on the galaxy 
activity. The CO cooling is maximum for the the '^CO(7-6) 
line and up for Henize 2-10 and the Cloverleaf QSO, while 
for IC 10, NGC 253, NGC 4038, M 83 and NGC 6946 
the turnover appears at the '-CO(6-5) line. For the center 
of the Milky Way and IC 342, the turnover is located at 
'■^CO(4-3) line. These lines are the main CO cooling lines, 
which shows the predominant role of the high-J CO lines 
(>'^CO(5-4)) in the molecular gas cooling. The shape of 
the CO cooling curve can be used as a diagnostic of the 
galaxy star forming activity. 

3 Whatever the morphological type, the C cooUng rates are 
lower than the CO cooling rates (the ratio varies from a 
factor of < 4.0 for IC 10 and IC 342, of > 10 for NGC 



253 and Henize 2-10 (reported in iBavet et al.ll2004l) . The 
factor are estimated to be ^ 20 for NGC 4038, Overlap, M 



7 While i-CO(3-2) and '^00(4-3) Unes are detected in 
star forming regions outside galaxy nuclei, we could not 
detect '-CO(6-5) emission there because of the lack of 
sensitivity. A similar analysis on spiral arm regions should 
be performed in the near future. 

8 For our galaxy nuclei sample, PDR solutions fitted to the 
CO lines, do not provide a good match to the fine structure 
line intensities of neutral carbon. 
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files either axisymmetric (See Eq. lA.lb or elliptical (See Eq. 
\A3 : 



1 JC^ 

fs(x,y) = fo — 7 exp 

no--; 



(A.l) 



fs{x,y) = fo- 



■ exp T- exp T- (A.2) 



where fo is a normalisation constant and where the Full 
Width Half Maximum (FWHM) is FWHM, = Icr, yfhH in the 
axisymmetric case, and FWHMj j - Icts^x Vln 2 et FWHM^ y 



— 2cr„ 



/In 2 in the elliptical case. The beam profile is also 



Gaussian, with a beam size FWHMof,., = 20-0^, Vln 2. The true 
source size (FWHMj) is obtained by performing a deconvolu- 
tion of the observed size using the known beam profile. 

The final line intensities is obtained by multiplying the ob- 
served signal by the following scaling factor, which depends on 
the source size, initial and final spatial resolution : 



fact — 



fi,. , +FWHM^ 

mh,imt s 

e^.f. +FWHM^ 



(A.3) 



for the axisymmetric case, and : 



fact 



fi, ^ , +FWHM2 

mb,mit 



9^, . , + FWHM?,. 



\ 0^ . ,. + FWHM? , \ 02 ,. + FWHM? , 

) mbjm -^^-^ y mhj in -^J 

for elliptical sources. 

Appendix B: Tables of data 



(A.4) 



Appendix A: Convolution formulae 

As CO lines have been observed with different spatial resolu- 
tions, the line intensities should be convolved to the same Unear 
resolution before performing meaningful comparisons. The fi- 
nal resolution of 21.9" has been chosen since it corresponds to 
the beaim size of the CSO when observing at 345 GHz (fre- 
quency of the '2CO(3-2)). To perform this convolution, we 
modelled emitting region of each galaxy using Gaussian pro- 
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Table B.l. Results from the gaussian fits derived from spectra of the Antennae (NGC 4038 et Overlap), Arp 220, Centaurus A, Henize 2-10, 
IC 10, IC 342, IRAS 10565, M 51, M 82, M 83, Mrk 231, NGC 253, NGC 891, NGC 1068, 3079, NGC 4736, NGC 6090 and NGC 6946. 
The observed lines ('^CO, '''CO and C) are indicated in the first column. For each source, the RA(1950) and DEC(1950) positions of the center 
(0",0") is listed in Table [T] The indication " corresponds to the spectra shown in this article (see Figs. [T]to I16t . 
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'-^€0(2-1) 


230.538 


30.50 


(-40,20) 


8.0±0.2 


1.0x10"'" 


2.5x10" 


18 






30.50 


(-40,10) 


12.3+0.2 


1.5x10"'° 


3.8x10" 


18 






30.50 


(-20,0) 


2.1+0.1 


2.6x10"" 


6.5x10" 


19 






30.50 


(0,-10) 


16.4+0.3 


2.1x10"'° 


5.1x10" 


18 






30.50 


(0,0)" 


22.3+0.2 


2.8x10"'° 


7.0x10" 


18 






30.50 


(0,10) 


19.1±0.2 


2.4x10"'° 


5.9x10" 


18 






30.50 


(10,0) 


26.1+0.2 


3.3x10"'° 


8.1x10" 


18 


...continued... 



...continued... 


Transition 


Freq 


beam 


offset 


j (T,„hdv) 


Intensity 

(Wm-^sr-i) 


Flux 






(GHz) 


size(") 


position(") 


(Kkms-i) 


(Wm-2 


) 






30.50 


(10,10) 


22.9±0.2 


2.9x10-'" 


7.1x10- 


i« 






30.50 


(10,20) 


11.4+0.2 


1.4x10-'" 


3.5x10- 


18 






30.50 


(20,0) 


13.0+0.3 


1.6x10-'° 


4.1x10- 


18 






30.50 


(30,10) 


4.5+0.3 


5.7x10-" 


1.4x10" 


18 


'"CO(3-2) 


345.796 


21.90 


(-45,20) 


3.0+0.8 


1.3x10-'° 


1.6x10- 


18 






21.90 


(-45,10) 


4.0+0.7 


1.7x10-'° 


2.2x10- 


18 






21.90 


(-35,10) 


5.9+0.7 


2.5x10-'° 


3.2x10" 


18 






21.90 


(-35,20) 


4.7+0.7 


2.0x10-'° 


2.5x10- 


18 






21.90 


(-5,0) 


11.5+0.9 


4.9x10'° 


6.2x10- 


18 






21.90 


(-5,10) 


10.0+0.9 


4.3x10"'° 


5.4x10- 


18 






21.90 


(5,-10) 


8.7+0.9 


3.7x10-'° 


4.7x10- 


18 






21.90 


(5,0)'' 


12.7+0.7 


5.3x10-'° 


6.8x10- 


18 






21.90 


(5,10) 


12.7+0.6 


5.4x10-'° 


6.9x10- 


18 






21.90 


(15,0) 


11.1+0.8 


4.7x10-'° 


6.0x10- 


18 






21.90 


(15,10) 


8.3+0.8 


3.5x10-'° 


4.5x10- 


18 


i^CO(3-2) 


330.588 


21.90 


(-10,0) 


0.8+0.4 


3.1x10-" 


4.0x10- 


19 






21.90 


(0,0)" 


0.9+0.1 


3.3x10-" 


4.3x10- 


19 






21.90 


(0,10) 


1.2+0.2 


4.6x10-" 


5.8x10- 


19 






21.90 


(0,20) 


0.6+0.3 


2.1x10-" 


2.7x10- 


19 






21.90 


(0,30) 


2.8+0.6 


1.0x10'° 


1.3x10- 


18 






21.90 


(10,0) 


0.9+0.3 


3.3x10-" 


4.3x10- 


19 


i-CO(4-3) 


461.041 


14.55 


(-40,10) 


4.1+0.6 


4.1x10-'° 


2.3x10- 


18 






14.55 


(-15,10) 


5.6+1.9 


5.6x10-'° 


3.2x10- 


18 






14.55 


(-15,20) 


3.9+2.3 


3.9x10-'° 


2.2x10- 


18 






14.55 


(-5,10) 


6.3+0.9 


6.3x10-'° 


3.5x10- 


18 






14.55 


(5,-10) 


8.8+1.9 


8.8x10-'° 


5.0x10- 


18 






14.55 


(5,0)" 


10.4+1.2 


LOxlO-** 


5.9x10- 


18 






14.55 


(5,10) 


11.3+1.0 


LlxlO-** 


6.4x10- 


18 






14.55 


(5,20) 


7.0+1.2 


7.0x10-'° 


4.0x10- 


18 






14.55 


(15,0) 


11.6+1.2 


1.2x10-** 


6.6x10- 


18 






14.55 


(15,10) 


10.7+1.9 


1.1x10-** 


6.1x10- 


18 






14.55 


(15,20) 


6.9+2.0 


6.9x10-'° 


3.9x10- 


18 


i2C0(6-5) 


691.473 


10.60 


(0,30) 


4.3+1.3 


1.5x10-'' 


4.4x10- 


18 






10.60 


(5,0) 


3.1+5.0 


l.OxlO-** 


3.1x10- 


18 






10.60 


(5,10)" 


11.2+1.7 


3.8x10-** 


1.1x10- 


17 


i2C0(7-6) 


806.652 


8.95 


(0,30)" 


3.6+0.8 


2.0x10-^ 


4.1x10- 


18 


IC 342 














CI(^P2-^Pl) 


809.902 


8.95 


(-5,-5) 


29.8+3.9 


1.6x10-** 


3.5x10- 


1/ 






8.95 


(-5,0) 


35.3+4.6 


1.9x10-** 


4.1x10- 


17 






8.95 


(0,0)" 


35.1+3.7 


1.9x10-^ 


4.1x10- 


17 






8.95 


(0,5) 


25.8+3.7 


1.4x10-** 


3.0x10- 


17 






8.95 


(5,-5) 


14.9+4.6 


8.1x10-** 


1.7x10" 


17 






8.95 


(5,0) 


49.8+6.3 


2.7x10-^ 


5.8x10" 


17 






8.95 


(5,5) 


31.8+6.6 


1.7x10-** 


3.7x10- 


17 






8.95 


(10,10) 


9.0+2.2 


4.9x10-** 


1.0x10- 


17 


'-CO(3-2) 


345.796 


21.90 


(-10,-10) 


72.9+6.4 


3.1x10-** 


3.9x10- 


17 






21.90 


(-10,0) 


70.1+3.2 


3.0x10-** 


3.8x10- 


17 






21.90 


(-10,10) 


51.0+5.1 


2.2x10-** 


2.8x10- 


17 






21.90 


(0,-10) 


56.1+7.1 


2.4x10-** 


3.0x10- 


17 






21.90 


(0,0)" 


109.8+2.5 


4.6x10-^ 


5.9x10- 


17 






21.90 


(0,10) 


97.0+3.5 


4.1x10-** 


5.2x10- 


17 






21.90 


(0,20) 


65.5+4.4 


2.8x10-** 


3.5x10" 


17 


...continued... 



...continued... 


Transition 


Freq 


beam 


offset 


j (T,„hdv) 


Intensity 

(Wm-^sr-i) 


Flux 






(GHz) 


size(") 


position(") 


(Kkms-i) 


(Wm-2) 






21.90 


(10,0) 


95.3±7.7 


4.0x10"'' 


5.1x10- 


1/ 






21.90 


(10,10) 


69.2+7.4 


2.9x10-'* 


3.7x10- 


17 






21.90 


(10,20) 


68.0+3.2 


2.9x10-'' 


3.7x10- 


17 


'^CO(3-2) 


330.588 


21.90 


(-10,0) 


6.7+0.5 


2.5x10-'° 


3.2x10- 


18 






21.90 


(0,-10) 


3.1+0.3 


1.1x10-'° 


1.4x10- 


18 






21.90 


(0,0)" 


19.9+2.2 


7.4x10-'° 


9.4x10- 


18 






21.90 


(0,10) 


9.0+0.4 


3.3x10-'° 


4.3x10- 


18 






21.90 


(10,0) 


2.9±0.7 


1.1x10-'° 


1.4x10- 


18 


'2cO(4-3) 


461.041 


14.55 


(0,0) 


110.8+12.6 


1.1x10-** 


6.3x10- 


17 


'2cO(6-5) 


691.473 


10.60 


(-5,-10) 


19.6+16.6 


6.6x10-'' 


2.0x10- 


17 






10.60 


(-5,-5) 


84.1+9.5 


2.8x10-** 


8.5x10- 


17 






10.60 


(-5,0) 


72.9+13.1 


2.5x10-^ 


7.4x10- 


17 






10.60 


(-5,5) 


23.7+6.4 


8.0x10-'' 


2.4x10- 


17 






10.60 


(0,0)" 


54.6+3.3 


1.8x10-^ 


5.5x10- 


17 






10.60 


(0,5) 


73.3+5.3 


2.5x10-** 


7.4x10- 


17 






10.60 


(0,10) 


54.8±9.6 


1.9x10-** 


5.5x10- 


17 






10.60 


(5,0) 


108.1+7.0 


3.7x10-** 


1.1x10- 


16 






10.60 


(5,5) 


90.9+5.5 


3.1x10-** 


9.2x10- 


17 






10.60 


(5,10) 


80.9+9.6 


2.7x10-** 


8.2x10- 


17 






10.60 


(5,15) 


17.5±4.3 


5.9x10-'' 


1.8x10- 


17 






10.60 


(10,5) 


57.0+7.1 


1.9x10-^ 


5.8x10- 


17 






10.60 


(10,10) 


81.6±7.3 


2.8x10-^ 


8.3x10- 


17 


i2C0(7-6) 


806.652 


8.95 


(-5,-5) 


39.8+3.4 


2.1x10-'^ 


4.6x10- 


17 






8.95 


(-5,0) 


33.0+6.9 


1.8x10-** 


3.8x10- 


17 






8.95 


(0,-5) 


25.1+5.9 


1.3x10-^ 


2.9x10- 


17 






8.95 


(0,0)" 


30.2+2.8 


1.6x10-** 


3.5x10- 


17 






8.95 


(0,5) 


33.3+4.1 


1.8x10-** 


3.8x10- 


17 






8.95 


(5,-5) 


19.0+5.6 


1.0x10-** 


2.2x10- 


17 






8.95 


(5,0) 


40.3±3.5 


2.2x10-** 


4.6x10- 


17 






8.95 


(5,5) 


29.7+4.5 


1.6x10-** 


3.4x10- 


17 






8.95 


(10,5) 


17.7+3.9 


9.5x10-'' 


2.0x10- 


17 


IRAS 10565 














CI(^Pi-^Po) 


492.162 


14.55 


(0,0)" 


11.4±2.1 


1.4x10-'' 


7.8x10- 


18 


'^CO(3-2) 


345.796 


21.90 


(0,0)" 


12.9±2.1 


5.5x10-'° 


7.0x10- 


18 


M51 














CiePi-'Po) 


492.162 


14.55 


(0,0)" 


15.1±2.6 


1.8x10-'' 


1.0x10- 


1/ 






14.55 


(-24,-24)" 


12.8+1.4 


1.6x10-^ 


8.8x10- 


18 






14.55 


(0,12)" 


8.1 + 1.5 


9.9x10-'° 


5.6x10- 


18 


'-CO(4-3) 


461.041 


14.55 


(-84,-84) 


5.5+1.1 


5.5x10-'° 


3.1x10- 


18 






14.55 


(-24,-24)" 


32.0+3.0 


3.2x10-'' 


1.8x10- 


17 






14.55 


(0,0)" 


36.7+2.4 


3.7x10-'' 


2.1x10- 


17 






14.55 


(0,12)" 


16.0+1.8 


1.6x10-^ 


9.1x10- 


18 


M82 














'-CO(2-l) 


230.538 


30.50 


(-10,-10) 


131.8+4.2 


1.7x10-^ 


4.1x10- 


17 






30.50 


(-10,0) 


267.7+3.5 


3.4x10-'' 


8.3x10- 


17 






30.50 


(-10,10) 


352.9+4.1 


4.4x10-'' 


1.1x10- 


16 






30.50 


(0,-10) 


148.9+4.0 


1.9x10-^ 


4.6x10- 


17 






30.50 


(0,0)" 


286.4+2.2 


3.6x10-'' 


8.9x10- 


17 






30.50 


(0,10) 


394.0+5.0 


4.9x10-'^ 


1.2x10- 


16 






30.50 


(10,-10) 


131.2+4.4 


1.6x10-'' 


4.1x10- 


17 






30.50 


(10,0) 


293.8+4.3 


3.7x10-'^ 


9.1x10- 


17 


...continued... 



...continued... 


Transition 


Freq 


beam 


offset 


j (T,„bdv) 


Intensity 

(Wm-^sr-i) 


Flux 






(GHz) 


size(") 


position(") 


(Kkms-i) 


(Wm-2) 






30.50 


(10,10) 


408.2±5.4 


5.1x10"'' 


1.3x10- 


16 






30.50 


(10,20) 


389.7±4.5 


4.9x10-'' 


1.2x10- 


16 






30.50 


(20,20) 


328.0+4.0 


4.1x10-^ 


1.0x10- 


16 






30.50 


(30,20) 


238.9±3.6 


3.0x10-'' 


7.4x10- 


17 


'"CO(6-5) 


691.473 


10.60 


(-10,-10) 


218.7±16.4 


7.4x10-** 


2.2x10- 


16 






10.60 


(-10,0) 


154.4±12.7 


5.2x10-** 


1.6x10- 


16 






10.60 


(-10,10) 


111.1±23.9 


3.8x10-** 


1.1x10- 


16 






10.60 


(-5,-5) 


170.0+28.2 


5.7x10-** 


1.7x10- 


16 






10.60 


(-5,0) 


202.6+19.9 


6.9x10-** 


2.0x10- 


16 






10.60 


(-5,5) 


252.5+25.1 


8.5x10-** 


2.6x10- 


16 






10.60 


(0,-10) 


194.8+30.2 


6.6x10-^ 


2.0x10- 


16 






10.60 


(0,-5) 


346.8+21.4 


1.2x10-^ 


3.5x10- 


16 






10.60 


(0,0)" 


263.0+21.0 


8.9x10-** 


2.7x10- 


16 






10.60 


(0,5) 


212.0+20.6 


7.2x10-^ 


2.1x10- 


16 






10.60 


(0,10) 


90.8+23.7 


3.1x10-^ 


9.2x10- 


17 






10.60 


(5,-5) 


188.2+26.6 


6.4x10-^ 


1.9x10- 


16 






10.60 


(5,0) 


280.0+19.6 


9.5x10-^ 


2.8x10- 


16 






10.60 


(5,5) 


366.2+28.2 


1.2x10-^ 


3.7x10- 


16 






10.60 


(10,-10) 


211.9+15.2 


7.2x10-** 


2.1x10- 


16 






10.60 


(10,0) 


282.1 + 19.2 


9.5x10-** 


2.9x10- 


16 






10.60 


(10,10) 


219.7+15.8 


7.4x10-** 


2.2x10- 


16 


M83 














'-^€0(2-1) 


230.538 


30.50 


(-10,0) 


62.4+1.3 


7.8x10-'" 


1.9x10- 


1/ 






30.50 


(0,0)" 


67.4+2.2 


8.4x10-'° 


2.1x10- 


17 






30.50 


(0,10) 


65.5+1.4 


8.2x10-"^ 


2.3x10- 


17 






30.50 


(0,-10) 


101.8+1.5 


1.3x10-'' 


3.2x10- 


17 






30.50 


(10,0) 


72.4+1.4 


9.1x10-"^ 


2.2x10- 


17 


i-CO(3-2) 


345.796 


21.90 


(-90,-60)" 


21.5+1.4 


9.1x10-'° 


1.2x10- 


17 






21.90 


(0,0)" 


99.9+1.6 


4.2x10-'' 


5.4x10- 


17 






21.90 


(90,30)" 


17.9+1.4 


7.6x10-'° 


9.7x10- 


18 


i^^CO(3-2) 


330.588 


21.90 


(-90,-60)" 


1.5+0.5 


5.5x10-" 


7.0x10" 


19 






21.90 


(0,0)" 


11.8+1.3 


4.3x10-'° 


5.5x10- 


18 






21.90 


(0,10) 


11.2+1.2 


4.3x10-'° 


5.3x10- 


18 






21.90 


(90,30)" 


1.8+0.3 


6.8x10" 


8.7x10- 


19 


i2C0(4-3) 


461.041 


14.55 


(-90,-60)" 


14.8+3.0 


1.5x10-'' 


8.4x10- 


18 






14.55 


(-5,-5) 


150.1+21.9 


1.5x10-^ 


8.5x10- 


17 






14.55 


(-5,5) 


136.1 + 13.2 


1.4x10-** 


7.7x10- 


17 






14.55 


(0,0)" 


159.0+16.4 


1.6x10-** 


9.0x10- 


17 






14.55 


(0,15) 


82.5+8.5 


8.3x10-'' 


4.7x10- 


17 






14.55 


(5,-5) 


114.1+9.4 


1.1x10-** 


6.4x10- 


17 






14.55 


(5,5) 


286.2+22.7 


2.9x10-** 


1.6x10- 


16 






14.55 


(10,10) 


112.3 + 13.3 


1.1x10-^ 


6.3x10- 


17 






14.55 


(15,0) 


82.4+16.5 


8.3x10-'' 


4.7x10- 


17 






14.55 


(100,30)" 


19.9+2.7 


2.0x10-'' 


1.1x10- 


17 


i2C0(6-5) 


691.473 


10.60 


(-10,0) 


134.8 + 13.7 


4.6x10-'^ 


1.4x10- 


16 






10.60 


(-10,5) 


65.6+15.5 


2.2x10-** 


6.6x10- 


17 






10.60 


(-5,-5) 


98.6+11.9 


3.3x10-^ 


1.0x10- 


17 






10.60 


(-5,0) 


317.5+20.6 


1.1x10-^ 


3.2x10- 


16 






10.60 


(-5,5) 


225.8+18.8 


7.6x10-** 


2.2x10- 


16 






10.60 


(-5,10) 


251.7+25.0 


8.5x10-** 


2.5x10- 


16 






10.60 


(0,-10) 


171.7+11.8 


5.8x10-^ 


1.7x10- 


16 






10.60 


(0,-5) 


162.0+12.8 


5.5x10-^ 


1.6x10- 


16 


...continued... 



...continued... 


Transition 


Freq 


beam 


offset 


j (T,„hdv) 


Intensity 

(Wm-^sr-i) 


Flux 






(GHz) 


size(") 


position(") 


(Kkms-i) 


(Wm"2 


) 






10.60 


(0,0)" 


235.3±14.5 


8.0x10"' 


2.4x10" 


16 






10.60 


(0,5) 


217.7±20.6 


7.4x10-** 


2.2x10" 


16 






10.60 


(0,10) 


135.2+14.6 


4.6x10"'^ 


1.4x10" 


16 






10.60 


(5,0) 


81.0+11.3 


2.7x10-^ 


8.2x10" 


17 






10.60 


(5,5) 


224.5+21.3 


7.6x10"** 


2.3x10" 


16 






10.60 


(5,10) 


124.0+15.8 


4.2x10"** 


1.3x10" 


16 


MRK 231 














CI(^Pi-^Po) 


492.162 


14.55 


(0,0)" 


7.5+1.3 


9.2x10"'" 


5.2x10" 


IX 


'-^€0(2-1) 


230.538 


30.50 


(0,0)" 


11.9+0.5 


1.5x10"'" 


3.7x10" 


1« 


i2C0(3-2) 


345.796 


21.90 


(0,0)" 


8.7+0.7 


3.7x10"'° 


4.7x10" 


18 


NGC 891 














CI(^Pi-^Po) 


492.162 


14.55 


(0,0) 


11.5+1.0 


1.4x10"'' 


7.9x10" 


18 






14.55 


(57,1) 


13.3+1.3 


1.6x10"^ 


9.2x10" 


18 






14.55 


(80,-6) 


14.9+1.7 


1.8x10"^ 


1.0x10" 


17 






14.55 


(120,-6) 


10.8+0.9 


1.3x10"'' 


7.4x10" 


18 


'^CO(2-l) 


220.39 


30.50 


(0,0) 


4.6+0.3 


5.8x10"" 


1.4x10" 


18 






30.50 


(9.3,14.5) 


5.3+0.4 


6.7x10"" 


1.7x10" 


18 






30.50 


(17.1,32.9) 


3.9+0.7 


4.9x10"" 


1.2x10" 


18 






30.50 


(25,51) 


2.7+0.2 


3.4x10"" 


8.4x10" 


19 






30.50 


(32.8,69.7) 


2.7+0.4 


3.4x10"" 


8.5x10" 


19 






30.50 


(40.6,88.1) 


2.1+0.2 


2.6x10"" 


6.4x10" 


19 






30.50 


(48.4,107) 


4.9+0.6 


6.2x10"" 


1.5x10" 


18 






30.50 


(56.2,125) 


2.4+0.2 


3.0x10"" 


7.3x10" 


19 






30.50 


(66.3,164) 


4.2+0.5 


5.2x10"" 


1.3x10" 


19 






30.50 


(82,201) 


1.4+0.2 


1.7x10"" 


4.2x10" 


19 


i-CO(3-2) 


345.796 


21.90 


(0,0) 


11.0+0.3 


4.7x10"'" 


6.0x10" 


18 






21.90 


(7.8,18.4) 


11.1+0.6 


4.7x10'" 


6.0x10" 


18 






21.90 


(9.3,14.5) 


10.2+0.3 


4.3x10"'" 


5.5x10" 


18 






21.90 


(25,51) 


10.1+0.4 


4.3x10"'" 


5.5x10" 


18 






21.90 


(32.8,69.7) 


14.0+0.5 


5.9x10"'" 


7.6x10" 


18 






21.90 


(40.6,88.1) 


13.2+0.4 


5.6x10"'" 


7.1x10" 


18 






21.90 


(48.4,106.5) 


13.8+2.0 


5.8x10"'" 


7.5x10" 


18 






21.90 


(56.2,125) 


3.2+0.3 


1.3x10"'" 


1.7x10" 


18 


NGC 3079 














CiePi-'Po) 


492.162 


14.55 


(0,0) 


59.6+3.8 


7.3x10"'' 


4.1x10" 


17 






14.55 


(0,0) 


56.7+3.2 


6.9x10"^ 


3.9x10" 


17 


'^CO(3-2) 


345.796 


21.90 


(0,0)" 


88.7+2.2 


3.7x10"'' 


4.8x10" 


17 


NGC 4736 














CI(^Pi-^Po) 


492.162 


14.55 


(0,0)" 


8.3+2.0 


1.0x10"'' 


5.7x10" 


18 


'^CO(2-l) 


220.399 


30.50 


(0,0)" 


2.2+0.6 


2.4x10"" 


5.9x10" 


19 






30.50 


(40,0) 


0.9+0.2 


1.1x10"" 


2.5x10" 


19 


i2C0(3-2) 


345.796 


21.90 


(0,0)" 


32.0+2.3 


1.4x10"'' 


1.7x10" 


17 






21.90 


(40,0) 


14.8+0.8 


6.3x10"'" 


8.0x10" 


18 


i3CO(3-2) 


330.588 


21.90 


(0,0)" 


0.9+0.4 


3.2x10"" 


4.1x10" 


19 


NGC 6090 














CI(^Pi-^Po) 


492.162 


14.55 


(0,0)" 


5.5+0.9 


6.8x10"'" 


3.8x10" 


18 


'-CO(2-l) 


230.538 


30.50 


(0,0)" 


13.8+1.1 


1.7x10"'" 


4.3x10" 


18 


'-CO(3-2) 


345.796 


21.90 


(0,0)" 


38.6+1.6 


1.6x10"^ 


2.1x10" 


17 


NGC 6946 














CI(^Pi-^Po) 


492.162 


14.55 


(0,0)" 


20.0+2.4 


2.4x10"'' 


1.4x10" 


17 






14.55 


(20,0) 


24.6+3.0 


3.0x10"^ 


1.7x10" 


17 
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Transition 


Freq 


beam 


offset 


j (T,„hdv) 


Intensity 

(Wm-^sr-i) 


Flux 




(GHz) 


size(") 


position(") 


(Kkms-i) 


(Wm-2) 






14.55 


(60,0) 


9.2±2.2 


1.1x10"'' 


6.4x10-'** 






14.55 


(110,100)" 


4.3+0.9 


5.3x10-'° 


3.0x10-'** 






14.55 


(150,-20)" 


4.0±0.9 


4.8x10-'° 


2.7x10-'** 


'-CO(2-l) 


230.538 


30.50 


(0,0)" 


104.7+2.8 


1.3x10"' 


3.3x10-" 






30.50 


(90,90) 


6.9+2.0 


8.7x10-" 


2.2x10-'** 






30.50 


(90,120) 


6.8+1.3 


8.5x10" 


2.1x10-'^ 






30.50 


(90,150) 


5.3±1.5 


6.6x10" 


1.6x10-'^ 






30.50 


(120,-30) 


10.8+2.8 


1.4x10-'° 


3.3x10-'* 






30.50 


(120,0) 


4.8+1.3 


6.0x10" 


1.5x10-'* 






30.50 


(150,-30) 


8.3+1.7 


1.0x10-'° 


2.6x10"'* 






30.50 


(150,0) 


6.5+1.4 


8.2x10"" 


2.0x10-'* 






30.50 


(180,-30) 


5.3+1.8 


6.7x10-" 


1.7x10"'* 


i2C0(3-2) 


345.796 


21.90 


(0,-10) 


90.0+2.0 


3.8x10-'^ 


4.9x10"''' 






21.90 


(0,0)" 


113.3+2.8 


4.7x10-'^ 


6.0x10-'^ 






21.90 


(10,0) 


134.9+2.1 


5.7x10-'^ 


7.3x10-" 






21.90 


(100,90) 


7.2+1.3 


3.1x10-'° 


3.9x10"'* 






21.90 


(100,100) 


9.7+2.1 


4.1x10-'° 


5.3x10-'* 






21.90 


(100,110)" 


6.6+1.3 


2.8x10-'° 


3.6x10-'* 






21.90 


(110,80) 


2.9+0.9 


1.2x10-'° 


1.5x10"'* 






21.90 


(110,100) 


14.7+1.7 


6.2x10-'° 


7.9x10"'* 






21.90 


(110,110) 


19.0+2.4 


8.0x10-'° 


1.0x10"''' 






21.90 


(140,-20) 


33.7+3.4 


1.4x10-'* 


1.8x10"''' 






21.90 


(150,-30) 


45.11 + 1.7 


1.9x10-'' 


2.4x10-''' 






21.90 


(150,-20)" 


9.0+1.6 


3.8x10-'° 


4.9x10-'* 






21.90 


(150,-10) 


10.2+1.7 


4.3x10-'° 


5.5x10"'* 






21.90 


(150,0) 


8.1 + 1.1 


3.4x10-'° 


4.4x10-'* 






21.90 


(160,-20) 


9.2+1.5 


3.9x10-'° 


5.0x10-'* 


i3CO(3-2) 


330.588 


21.90 


(0,0)" 


7.9+0.8 


2.9x10-'° 


3.8x10-'* 


i-CO(4-3) 


461.041 


14.55 


(0,0)" 


194.1+5.7 


1.9x10-** 


LlxlO-"" 






14.55 


(110,100)" 


9.2+1.3 


9.2x10-'° 


5.2x10-'* 






14.55 


(150,-20)" 


14.4+2.3 


1.4x10-'' 


8.1x10-'* 


i2C0(6-5) 


691.473 


10.60 


(0,0)" 


70.1+10.4 


2.4x10-^ 


7.1x10-" 






10.60 


(0,10) 


51.8+12.7 


1.8x10-** 


5.2x10-" 






10.60 


(5,-10) 


111.4+23.0 


3.8x10-^ 


LlxlO-'*- 






10.60 


(5,0) 


90.7+9.8 


3.1x10-^ 


9.2x10-''' 






10.60 


(5,5) 


102.4+16.8 


3.5x10-** 


1.0x10-"' 






10.60 


(5,10) 


83.3+11.7 


2.8x10-** 


8.4x10-" 






10.60 


(5,15) 


129.4+55.1 


4.4x10-** 


1.3x10-"' 






10.60 


(10,0) 


152.4+23.7 


5.2x10-** 


1.5x10'*' 






10.60 


(10,10) 


181.5+53.7 


6.1x10-** 


1.8x10"'*' 



END 
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Table B.2. Results of the gaussian fits derived from spectra of the Antennae (NGC 4038 et Overlap), Henize 2-10, IC 10, IC 342, M 83, NGC 
253 and NGC 6946. See caption of Table IB. II Data listed in this table are either from literature or from our sample of observations. 



Transition 


Freq 


beam 


offset 


/ (Tmhdv) 


Intensity 


Flux 




Ref. 




(GHz) 


size(") 


position(") 


(Kkms"') 


(Wm^^sr') 


(Wm-2) 




Antennae 


















NGC4038 


















CI(^Pi-^Po) 


492.162 


14.55 


(0,0) 


13.7+2.1 


1.7x10"^ 


9.4x10- 


18 


1 






21.90 


(0,0) 


7.6±1.2 


9.2x10-'" 


1.2x10- 


17 


1* 


i^CO(l-O) 


115.271 


43.00 


(0,0) 


21.0±2.0 


3.3x10-" 


1.6x10- 


18 


2 






6.75 


(0,0) 


415.0±83.0 


6.5x10-'" 


7.9x10- 


19 


3 






21.90 


(0,0) 


112.1+22.4 


1.8x10-1" 


2.2x10- 


18 


3 






55.00 


(0,0) 


19.5±3.9 


3.1x10-" 


2.5x10- 


18 


6 






15.00 


(0,0) 


80.0+16.0 


1.3x10-1" 


7.5x10- 


19 


4 






21.90 


(0,0) 


45.9+9.2 


7.2x10-11 


9.2x10- 


19 


4* 


i2CO(2-l) 


230.538 


30.50 


(0,0) 


35.6+1.3 


4.5x10-'" 


1.1x10- 


17 


5 






20.00 


(0,0) 


59.1 + 11.8 


7.4x10-'" 


7.9x10- 


18 


4 






21.90 


(0,0) 


51.3+10.3 


6.4x10-'" 


8.2x10- 


18 


4* 






30.50 


(0,0) 


31.3+6.3 


3.9x10-'" 


9.7x10- 


18 


7 


'2CO(3-2) 


345.796 


21.90 


(0,0) 


37.2+1.7 


1.6x10-'' 


2.0x10- 


17 


5* 






14.00 


(0,0) 


56.1 + 11.2 


2.4x10-'' 


1.2x10- 


17 


4 






21.90 


(0,0) 


29.4+5.9 


1.2x10-'' 


1.6x10- 


17 


4 


'2C0(4-3) 


461.041 


14.55 


(0,0) 


45.3+3.5 


4.5x10-'^ 


2.6x10- 


17 


5 






21.90 


(0,0) 


25.0+1.9 


2.5x10-'^ 


3.2x10- 


7 


5* 


i2C0(6-5) 


691.473 


10.60 


(3,3) 


85.5+10.4 


2.9x10-** 


8.7x10- 


17 


5 






21.90 


(3,3) 


32.8+4.0 


1.1x10-** 


1.4x10- 


6 


5* 


'^CO(7-6) 


806.652 


8.95 


(3,3) 


4.9+1.3 


2.6x10-'' 


5.6x10- 


8 


5 






21.90 


(3,3) 


1.6+0.4 


8.7x10 


1.1x10- 


7 


5* 


^■^CO(2-l) 


220.399 


20.00 


(0,0) 


2.4+0.5 


2.7x10-" 


2.8x10- 


19 


4* 






21.90 


(0,0) 


2.1+0.4 


2.3x10-" 


2.9x10- 


19 


4 






30.50 


(0,0) 


2.4+0.5 


2.6x10-" 


6.5x10- 


19 


7 


'^CO(3-2) 


330.588 


21.90 


(0,0) 


2.1 + 1.2 


7.6x10-" 


9.8x10- 


19 


5* 






14.00 


(0,0) 


4.5+0.9 


1.7x10-'" 


8.7x10- 


19 


4 






21.90 


(0,0) 


2.4+0.5 


8.7x10-" 


1.1x10- 


18 


4 


Antennae 


















Overlap 


















CI(^Pi-^Po) 


492.162 


14.55 


(0,0) 


11.2+1.5 


1.4x10-'' 


7.7x10- 


18 


1 






21.90 


(0,0) 


5.8+0.8 


7.1x10-'" 


9.1x10- 


18 


1* 


'-CO(l-O) 


115.271 


6.75 


(0,0) 


721.4+144.3 


1.1x10-'^ 


1.4x10- 


18 


2 






21.9 


(0,0) 


163.4+32.7 


2.6x10-'" 


3.3x10- 


18 


2 






43.00 


(0,0) 


22.0+2.0 


3.4x10-" 


1.7x10- 


18 


3 






15.00 


(0,0) 


128.0+25.6 


2.0x10-1" 


1.2x10- 


18 


4 






21.90 


(0,0) 


70.0+14.0 


1.1x10-'" 


1.4x10- 


18 


4* 






55.00 


(0,0) 


19.2+3.8 


3.0x10-" 


2.4x10- 


18 


6 


i2C0(2-l) 


230.538 


30.50 


(0,0) 


58.3+2.6 


7.3x10-'" 


1.8x10- 


17 


5 






20.00 


(0,0) 


87.9+17.6 


1.1x10-'^ 


1.2x10- 


17 


4 






21.90 


(0,0) 


75.4+15.1 


9.5x10-'" 


1.2x10- 


17 


4* 






30.50 


(0,0) 


61.2+12.2 


7.7x10-'" 


1.9x10- 


17 


7 


i2CO(3-2) 


345.796 


21.90 


(0,0) 


56.8+3.5 


2.4x10-'^ 


3.1x10- 


17 


5* 






14.00 


(0,0) 


70.3+14.1 


3.0x10-'' 


1.6x10- 


17 


4 






21.90 


(0,0) 


34.8+7.0 


1.5x10-'' 


1.9x10- 


17 


4 


i2CO(4-3) 


461.041 


14.55 


(0,0) 


45.1+2.5 


4.5x10-'' 


2.5x10- 


17 


5 






21.90 


(0,0) 


23.6+1.3 


2.4x10-'' 


3.0x10- 


17 


5* 


i2C0(6-5) 


691.473 


10.60 


(3,3) 


97.0+37.1 


3.3x10-** 


9.8x10- 


17 


5 
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Transition 


Freq 


beam 


offset 


j {T„,hdv) 


Intensity 


Flux 


Ref. 




(GHz) 


size(") 


position(") 


(Kkms-i) 


(Wm-2sr-i) 


(Wm-2) 


a 






21.90 


(3,3) 


33.6±12.8 


1.1x10"' 


1.4x10-'" 


5* 


'2C0(7-6) 


806.652 


8.95 


(3,3) 


13.4+1.6 


7.2x10-'* 


1.5x10-'^ 


5 






21.90 


(3,3) 


3.9±0.5 


2.1x10-'' 


2.7x10-'^ 


5* 


'"CO(2-l) 


220.399 


20.00 


(0,0) 


5.4±1.1 


5.9x10-" 


6.3x10-''' 


4 






21.90 


(0,0) 


4.6+0.9 


5.1x10-" 


6.5x10-''' 


4* 






30.50 


(0,0) 


3.3±0.7 


3.6x10" 


9.0x10-''' 


7 


i3C0(3-2) 


330.588 


14.00 


(0,0) 


3.8±0.8 


1.4x10-'" 


7.3x10-''' 


4 






21.90 


(0,0) 


1.9±0.4 


7.0x10-" 


8.9x10-''' 


4* 


HENIZE 
















2-10 
















CI(^Pi-^Po) 


492.162 


14.55 


(0,0) 


4.2±0.8 


5.1x10-'" 


2.9x10-"* 


4et5 






21.90 


(0,0) 


2.5+0.5 


3.0x10-'" 


3.8x10-"* 


4et5* 


'^CO(l-O) 


115.271 


40.0 


(0,0) 


10.0+0.8 


1.6x10-" 


6.7x10-''' 


1 






21.90 


(0,0) 


27.3+2.2 


4.3x10-" 


5.5x10-'^ 


1 






55.00 


(0,0) 


4.9+0.2 


7.7xl0-'2 


6.2x10-'^ 


2 






21.90 


(0,0) 


24.1±1.0 


3.8x10-" 


4.8x10-''' 


2* 


'^CO(2-l) 


230.538 


21.00 


(0,0) 


17.3±1.4 


2.2x10-'" 


2.5x10-"* 


1* 






27.00 


(0,0) 


6.8±0.8 


8.5x10-" 


1.7x10-"* 


2 






21.90 


(0,0) 


9.4+1.1 


1.2x10-'" 


1.5x10-"* 


2 


'2CO(3-2) 


345.796 


21.90 


(0,0) 


11.5±2.3 


4.9xl0~'° 


6.2x10'^ 


5* 






21.00 


(0,0) 


23.2+2.1 


9.8x10-'" 


1.2x10-'^ 


1 






22.00 


(0,0) 


16.6+0.6 


7.0x10-'" 


9.0x10-"* 


3 


'2CO(4-3) 


461.041 


14.55 


(0,0) 


18.6±4.2 


1.9x10-'' 


1.1x10-'^ 


5 






21.90 


(0,0) 


10.9+2.4 


1.1x10-'' 


1.4x10-'^ 


5* 


'2C0(6-5) 


691.473 


10.60 


(0,0) 


15.7+3.1 


5.3x10-'' 


1.6x10-'^ 


5 






21.90 


(0,0) 


6.8+1.3 


2.3x10-^ 


2.9x10-'^ 


5* 


i2CO(7-6) 


806.652 


8.95 


(0,0) 


15.2+3.0 


8.2x10-^ 


1.7x10-'^ 


5 






21.90 


(0,0) 


5.8+1.2 


3.1x10-^ 


4.0x10-'^ 


5* 


'^CO(l-O) 


110.201 


57.00 


(0,0) 


0.3±0.1 


4.1x10-'' 


3.6x10-^" 


2 






21.90 


(0,0) 


1.6±0.5 


2.2xl0-'2 


2.8x10-2" 


2* 






40.00 


(0,0) 


<0.5 


< 6.8x10-'^ 


< 2.9x10-2" 


1 






21.90 


(0,0) 


< 


< xl0-'2 


< xlO-2" 


1 


i^^CO(2-l) 


220.399 


21.00 


(0,0) 


0.9±0.2 


9.9x10-'^ 


1.2x10-''' 


1* 


i^^CO(3-2) 


330.588 


21.90 


(0,0) 


1.8+0.5 


6.7x10-" 


8.5x10-''' 


5 






14.00 


(0,0) 


2.3+0.6 


8.5x10-" 


4.4x10"'^ 


1 






21.90 


(0,0) 


1.3+0.3 


4.8x10-" 


6.1x10-'^ 


1* 


IC 10 
















CI(^Pi-^Po) 


492.162 


14.55 


(5,10) 


5.5±0.6 


6.8x10-'" 


3.8x10-"* 


1 






21.90 


(5,10) 


3.3+0.4 


4.1x10-'" 


5.2x10-"* 


1* 






10.80 


(0,-3) 


2.2+0.5 


2.7x10-'" 


8.4x10-''' 


2 






21.90 


(0,-3) 


1.0+0.2 


1.2x10-'" 


1.6x10-"* 


2* 


CI(3p2-3Pi) 


809.902 


8.95 


(0,30) 


3.0+0.7 


1.6x10-'' 


3.5x10-"* 


3 






21.90 


(0,30) 


1.2±0.3 


6.6x10-'" 


8.4x10-"* 


3* 


'^CO(l-O) 


115.271 


22.0 


(5,10) 


22.8±4.6 


3.6x10-" 


4.6x10-''' 


4* 


i2C0(2-l) 


230.538 


30.50 


(10,0) 


26.1+0.2 


3.3x10-'" 


8.1x10-"* 


3 






21.90 


(10,0) 


43.7+0.3 


5.5x10-'" 


7.0x10-"* 


3 






12.50 


(5,0) 


33.9+6.8 


4.2x10-'" 


1.8x10-"* 


4 






21.90 


(5,0) 


17.5+3.5 


2.2x10-'" 


2.8x10-"* 


4 






22.00 


(2.2,10) 


15.7±3.1 


2.0x10-'" 


2.5x10-"* 


5* 






30.50 


(0,30) 


5.5±1.1 


6.9x10-" 


1.7x10-"* 


6 






21.90 


(0,30) 


9.2+1.8 


1.2x10-'" 


1.5x10-"* 


6 



...continued... 



...continued... 



Transition 


Freq 


beam 


offset 


j {T„,hdv) 


Intensity 


Flux 




Ref. 




(GHz) 


size(") 


position(") 


(Kkms-i) 


(Wm-2sr-i) 


(Wm-2) 


a 


'''CO(3-2) 


345.796 


21.90 


(5,10) 


12.7±0.6 


5.4x10"'" 


6.9x10- 




3* 






22.00 


(2.2,10) 


16.3+3.3 


6.9x10-1° 


8.9x10- 


18 


5 






13.20 


(0,-3) 


14.9+0.6 


6.3x10-'° 


2.9x10- 


18 


2 






21.90 


(0,-3) 


8.1+0.3 


3.4x10-'° 


4.4x10- 


18 


2 


'^CO(4-3) 


461.041 


14.55 


(15,0) 


11.6+1.2 


1.2x10-'* 


6.6x10- 


18 


3 






21.90 


(15,0) 


7.0+0.7 


7.0x10-'° 


8.9x10- 


18 


3* 


i2C0(6-5) 


691.473 


10.60 


(5,10) 


11.2+1.7 


3.8x10-^ 


1.1x10- 


17 


3 






21.90 


(5,10) 


5.1+0.8 


1.7x10-'* 


2.2x10- 


17 


3* 


'^CO(7-6) 


806.652 


8.95 


(0,30) 


3.6+0.8 


1.9x10-'* 


4.2x10- 


18 


3 






21.90 


(0,30) 


1.5+0.3 


7.9x10-'° 


1.0x10- 


17 


3* 


'^'CO(2-l) 


220.399 


30.50 


(0,30) 


0.4+0.1 


4.3x10-'" 


1.1x10- 


1^) 


1 






21.90 


(0,30) 


0.7+0.2 


7.3x10-'" 


9.4x10- 


20 


1 






22.00 


(2.2,10) 


1.0+0.2 


1.1x10-" 


1.4x10- 


19 


5* 


'3CO(3-2) 


330.588 


15.00 


(2.2,10) 


4.4+0.9 


1.6x10-'° 


9.6x10- 


19 


5 






21.90 


(2.2,10) 


2.7+0.6 


1.0x10-'° 


1.3x10- 


18 


5 






13.20 


(0,-3) 


2.5+0.6 


9.2x10-" 


4.3x10- 


19 


2 






21.90 


(0,-3) 


1.4+0.3 


5.0x10-" 


6.4x10- 


19 


2 






21.90 


(0,10) 


1.2+0.2 


4.6x10-" 


5.9x10- 


19 


3* 


IC 342 
















CI(^Pi-^Po) 


492.162 


10.00 


(0,0) 


54.0+6.0 


6.6x10-'' 


1.8x10- 


1/ 


1 






21.90 


(0,0) 


12.0+1.3 


1.5x10-'^ 


1.9x10- 


17 


1* 


CI(3P2-3P1) 


809.902 


8.95 


(5,0) 


49.8+6.3 


2.7x10-'* 


5.8x10- 


17 


2 






21.90 


(5,0) 


9.1 + 1.1 


4.9x10-'' 


6.3x10- 


17 


2* 


'-=€0(1-0) 


115.271 


21.00 


(0,0) 


213.2+2.1 


3.3x10-'° 


3.9x10- 


18 


3* 


'2C0(2-1) 


230.538 


14.00 


(0,0) 


324.3+64.8 


4.1x10-'* 


2.1x10- 


17 


3 






21.90 


(0,0) 


135.9+27.2 


1.7x10-'* 


2.2x10" 


17 


3* 






21.00 


(0,0) 


172.0+19.0 


2.2x10-'* 


2.5x10- 


17 


1 


'2CO(3-2) 


345.796 


14.00 


(0,0) 


186.0+23.0 


7.9x10-'* 


4.1x10- 


17 


1 






21.90 


(0,0) 


78.0+9.6 


3.3x10-'' 


4.2x10- 


17 


1 






21.90 


(0,0) 


109.8+2.5 


4.6x10-'* 


5.9x10- 


17 


2* 


i2CO(4-3) 


461.041 


11.00 


(0,0) 


209.0+21.0 


2.1x10-** 


6.7x10- 


17 


1 






21.90 


(0,0) 


55.5+5.6 


5.6x10-^ 


7.1x10- 


17 


1 






14.55 


(0,0) 


110.8+12.6 


1.1x10-** 


6.3x10- 


17 


2 






21.90 


(0,0) 


50.0+5.7 


5.0x10-'* 


6.4x10- 


17 


2* 


i2CO(6-5) 


691.473 


10.60 


(5,0) 


108.1+7.0 


3.7x10-^ 


1.1x10- 


16 


2 






21.90 


(5,0) 


26.8+1.7 


9.1x10-^ 


1.2x10- 


16 


2* 


'2CO(7-6) 


806.652 


8.95 


(5,0) 


40.3+3.5 


2.2x10-^ 


4.6x10- 


17 


2 






21.90 


(5,0) 


7.3+0.6 


3.9x10-'* 


5.0x10- 


17 


2* 


'"CO(l-O) 


110.201 


21.00 


(0,0) 


18.9+0.5 


2.6x10-" 


3.0x10- 


19 


3* 


'3CO(2-l) 


220.399 


14.00 


(0,0) 


29.7+1.1 


3.3x10-'° 


1.7x10- 


18 


3 






21.90 


(0,0) 


12.4+0.5 


1.4x10-'° 


1.7x10- 


18 


3* 






21.00 


(0,0) 


24.0+3.0 


2.6x10-'° 


3.1x10- 


18 


1 


'3CO(3-2) 


330.588 


14.00 


(0,0) 


17.1+2.0 


6.3x10-'° 


3.3x10- 


18 


1 






21.90 


(0,0) 


7.2+0.8 


2.6x10'° 


3.4x10- 


18 


1 






21.90 


(0,0) 


19.9+2.2 


7.4x10-'° 


9.4x10- 


18 


2* 


M83 
















CI(^^Pi-^Po) 


492.162 


10.00 


(5,5) 


78.0+15.6 


9.5x10-^^ 


2.5x10- 


1/ 


1 






21.90 


(5,5) 


18.7+3.7 


2.3x10-'* 


2.9x10- 


17 


1* 






21.00 


(0,0) 


55.0+8.0 


6.7x10-'* 


7.9x10- 


17 


2 


'-=€0(1-0) 


115.271 


24.00 


(0,0) 


87.7+2.1 


1.3x10-'° 


2.1x10- 


18 


3* 






16.00 


(0,0) 


80.0+16.0 


1.3x10-'° 


8.5x10- 


19 


4 



...continued... 



...continued... 



Transition 


Freq 


beam 


offset 


j {T„,hdv) 


Intensity 


Flux 




Ref. 




(GHz) 


size(") 


position(") 


(Kkms-1) 


(Wm-2sr-i) 


(Wm"2) 


a 






21.90 


(0,0) 


44.2+8.8 


6.9x10"" 


8.8x10" 


ly 


4 


i2C0(2-l) 


230.538 


24.00 


(0,0) 


98.1+0.8 


1.2x10-^ 


1.9x10" 


17 


3* 






21.00 


(0,0) 


261.0+52.2 


3.3x10-'* 


3.8x10" 


17 


2 






30.50 


(0,-10) 


101.8+20.4 


1.3x10"'' 


3.2x10" 


17 


5 


'^CO(3-2) 


345.796 


21.90 


(0,0) 


99.9+1.6 


4.2x10"'' 


5.4x10" 


17 


5* 






14.00 


(5,5) 


290.0+58.0 


1.2x10"** 


6.4x10" 


17 


1 






21.90 


(5,5) 


125.3+25.1 


5.3x10"^ 


6.8x10" 


17 


1 






21.00 


(0,0) 


167.0+15.0 


7.1x10"'' 


8.3x10" 


17 


2 


i2C0(4-3) 


461.041 


14.55 


(5,5) 


196.9+39.4 


2.0x10"** 


1.1x10" 


16 


5 






21.90 


(5,5) 


91.3+18.3 


9.1x10"'' 


1.2x10" 


16 


5* 






21.00 


(0,0) 


122.0+15.0 


1.2x10"^ 


1.4x10" 


16 


2 






11.00 


(5,5) 


357.0+71.4 


3.6x10"** 


1.2x10" 


16 


1 






21.90 


(5,5) 


100.5+20.1 


1.0x10"^ 


1.3x10" 


16 


1 


'2CO(6-5) 


691.473 


10.60 


(-5,0) 


317.5+20.6 


1.1x10"^ 


3.2x10" 


16 


5 






21.90 


(-5,0) 


83.9+5.4 


2.8x10"^ 


3.6x10" 


16 


5* 


'"CO(2-l) 


220.399 


21.00 


(0,0) 


28.5+3.0 


3.1x10""' 


3.7x10" 


IH 


2* 


i3CO(3-2) 


330.588 


21.90 


(0,0) 


11.8+1.3 


1.3x10"'" 


1.7x10" 


18 


5* 






14.00 


(0,0) 


22.3+1.0 


8.2x10"'° 


4.3x10" 


18 


2 






21.90 


(0,0) 


9.6+0.4 


3.6x10"'" 


4.5x10" 


18 


2 


NGC 253 
















CI(^Pi-^Po) 


492.162 


22.00 


(0,0) 


290.0+45.0 


3.5x10"** 


4.6x10" 


16 


10* 






10.20 


(0,0) 


575.0+115.0 


7.0x10"** 


1.9x10" 


16 


5 






21.90 


(0,0) 


210.5+42.1 


2.6x10"** 


3.3x10" 


16 


5 






43.00 


(0,0) 


98.0+19.6 


1.2x10"^ 


5.9x10" 


16 


5 






23.00 


(0,0) 


320.0+64.0 


3.9x10"** 


5.5x10" 


16 


4 


CI(3P2-3P1) 


809.902 


8.95 


(0,0) 


188.5+37.7 


1.0x10"^ 


2.2x10' 


16 


11 






21.90 


(0,0) 


59.8+12.0 


3.3x10"^ 


4.2x10" 


16 


11* 


'''CO(l-O) 


115.271 


43.00 


(0,0) 


343.0+68.6 


5.4x10"'" 


2.6x10" 


17 


2 






23.00 


(0,0) 


920.0+82.8 


1.4x10"'' 


2.0x10" 


17 


6* 


'2C0(2-1) 


230.538 


23.00 


(0,0) 


1062.0+116.8 


1.3x10"^ 


1.9x10" 


16 


7* 






21.00 


(0,0) 


926.0+185.2 


1.2x10"^ 


1.4x10" 


16 


2 


^2CO(3-2) 


345.796 


21.90 


(0,0) 


815.6+163.1 


3.4x10"^ 


4.4x10" 


16 


11* 






23.00 


(0,0) 


998.0+139.7 


4.2x10"** 


5.9x10" 


16 


7 






23.00 


(0,0) 


1194.0+238.8 


5.0x10"** 


7.1x10" 


16 


2 






14.00 


(0,0) 


1200.0+240.0 


5.1x10"** 


2.6x10" 


16 


5 






21.90 


(0,0) 


642.3+128.5 


2.7x10"** 


3.5x10" 


16 


5 






22.00 


(0,0) 


680.0+60.0 


2.9x10"** 


3.7x10" 


16 


8 


'2CO(4-3) 


461.041 


15.00 


(0,0) 


507.0+101.4 


5.1x10"** 


3.0x10" 


16 


3 






21.90 


(0,0) 


296.8+59.4 


3.0x10"^ 


3.8x10" 


16 


3 






22.00 


(0,0) 


1019.0+120.0 


1.0x10"^ 


1.3x10" 


15 


9* 






10.40 


(0,0) 


2160.0+432.0 


2.2x10"^ 


6.2x10" 


16 


5 






21.90 


(0,0) 


808.4+161.7 


8.1x10"** 


1.0x10" 


15 


5 


'2C0(6-5) 


691.473 


10.60 


(0,0) 


1394.0+278.8 


4.7x10"^ 


1.4x10" 


15 


11 






21.90 


(0,0) 


533.1 + 106.6 


1.8x10"^ 


2.3x10" 


15 


11* 






8/30 


(0,0) 


861+258.3 


2.9x10"^ 


5.0x10" 


15 


1 


i2CO(7-6) 


806.652 


8.95 


(0,0) 


810.2+162.0 


4.3x10"^ 


9.3x10" 


16 


11 






21.90 


(0,0) 


257.2+51.4 


1.4x10"^ 


1.8x10" 


15 


11* 






11.5/60 


(0,0) 


1370+411 


7.3x10"^ 


2.6x10" 


15 


10 


''CO(l-O) 


110.201 


23.00 


(0,0) 


80.0+8.0 


1.1x10"'" 


1.5x10" 


18 


7* 


^3CO(2-l) 


220.399 


23.00 


(0,0) 


82.0+9.8 


9.0x10"'" 


1.3x10" 


17 


7* 






21.00 


(0,0) 


104.0+20.8 


1.1x10"'' 


1.3x10" 


17 


2 



...continued... 



...continued... 


Transition 


Freq 


beam 


offset 


j (T,„hdv) 


Intensity 


Flux 




Ref. 




(GHz) 


size(") 


position(") 


(Kkms-1) 


(Wm-2sr-i) 


(Wm-2) 


a 


'^'CO(3-2) 


330.588 


23.00 


(0,0) 


90.0±12.6 


3.3x10"'' 


4.7x10- 


/ 








23.00 


(0,0) 


210.0+42.0 


7.8x10-^ 


1.1x10- 


6 


2 


NGC 6946 
















CI(^Pi-^Po) 


492.162 


21.00 


(0,0) 


44.0±8.0 


5.4x10"' 


6.3x10- 


/ 


1 






14.55 


(20,0) 


24.6+3.0 


3.0x10-5 


1.7x10- 


7 


2 






21.90 


(20,0) 


11.3±1.4 


1.4x10-'* 


1.8x10- 


7 


2* 


'''CO(l-O) 


115.271 


21.00 


(0,0) 


198.1+39.6 


3.1x10-'" 


3.6x10- 


« 


3* 






5.65 


(0,0) 


945.0+189.0 


1.5x10-'' 


1.3x10- 


8 


4 






21.90 


(0,0) 


93.5+18.7 


1.5x10-'" 


1.9x10- 


8 


4 






23.00 


(0,0) 


189.8+38.0 


3.0x10-'" 


4.2x10- 


8 


5 






17.00 


(0,0) 


227.4+45.5 


3.6x10'° 


2.7x10- 


8 


6 






21.90 


(0,0) 


140.2+28.0 


2.2x10-'° 


2.8x10- 


8 


6 






23.00 


(0,0) 


169.7+33.9 


2.7x10-'° 


3.7x10- 


8 


7 


'2CO(2-l) 


230.538 


30.50 


(0,0) 


104.7+2.8 


1.3x10-'' 


3.3x10- 


7 


8 






21.90 


(0,0) 


199.7+5.3 


2.5x10-5 


3.2x10- 


7 


8* 






21.00 


(0,0) 


222.0+20.0 


2.8x10-'' 


3.3x10- 


7 


1 






14.00 


(0,0) 


170.0+34.0 


2.1x10-'' 


1.1x10- 


7 


9 






21.90 


(0,0) 


73.0+14.6 


9.1x10-'° 


1.2x10- 


7 


9 






14.00 


(0,0) 


178.3+35.7 


2.2x10-'' 


1.2x10- 


7 


7 






21.90 


(0,0) 


76.5+15.3 


9.6x10'° 


1.2x10- 


7 


7 


i2C0(3-2) 


345.796 


22.00 


(0,0) 


129.1+25.2 


5.5x10-'' 


7.0x10- 


7 


3 






21.00 


(0,0) 


145.0+15.0 


6.1x10'' 


7.2x10- 


7 


1 






21.90 


(0,0) 


200.0+40.0 


8.5x10'' 


1.1x10- 


6 


10 






21.90 


(10,0) 


134.9+2.1 


5.7x10-'' 


7.3x10- 


7 


8* 






22.00 


(0,0) 


46.0+2.0 


1.9x10-'' 


2.5x10- 


7 


11 


'2CO(4-3) 


461.041 


17.00 


(0,0) 


82.4+16.5 


8.3x10-5 


6.4x10- 


7 


12 






21.90 


(0,0) 


50.8+10.2 


5.1x10-5 


6.5x10- 


7 


12 






21.00 


(0,0) 


112.0+11.0 


1.1x10-^ 


1.3x10- 


6 


1 






17.00 


(0,0) 


155.4+31.1 


1.6x10-^ 


1.2x10- 


6 


3 






21.90 


(0,0) 


95.8+19.2 


9.6x10-5 


1.2x10- 


6 


3 






14.55 


(0,0) 


194.1+38.8 


1.9x10-^ 


1.1x10- 


6 


8 






21.90 


(0,0) 


89.4+17.9 


9.0x10-5 


1.1x10- 


6 


8* 


i-CO(6-5) 


691.473 


10.60 


(5,-10) 


111.4+23.0 


3.8x10-** 


1.1x10- 


6 


8 






21.90 


(5,-10) 


29.1+6.0 


9.8x10-5 


1.3x10- 


6 


8* 


^^CO(l-O) 


110.201 


21.00 


(0,0) 


17.9+3.6 


2.4x10-" 


2.9x10- 


9 


3* 


i3CO(2-l) 


220.399 


12.00 


(0,0) 


20.7+4.1 


2.3x10-'° 


8.7x10- 


9 


3 






21.90 


(0,0) 


6.7+1.3 


7.4x10-" 


9.4x10- 


9 


3* 






14.00 


(0,0) 


22.2+3.0 


2.4x10-'° 


1.3x10- 


8 


1 






21.90 


(0,0) 


9.5+1.3 


1.0x10-'° 


1.3x10- 


8 


1 


"CO(3-2) 


330.588 


22.00 


(0,0) 


5.6+1.1 


6.1x10-" 


7.9x10- 


9 


3 






14.00 


(0,0) 


11.4+2.0 


4.2x10-'° 


2.2x10- 


8 


1 






21.90 


(0,0) 


4.9+0.9 


1.8x10-'° 


2.3x10- 


8 


1 






21.90 


(0,0) 


28.0+5.6 


1.0x10-5 


1.3x10- 


7 


10 






21.90 


(0,0) 


7.9+0.8 


2.9x10-'° 


3.7x10- 


8 


8* 
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Refe ren ces: Pour NGC 4038 : 1: Gerin & Phillipsl (llOOOl) but spec tra have been an aly zed again; 2 



Aalto et al. 



( Il995l) ; 3: [Stanford et al.'('l990"); 4:'Zhu et al.' ('2003"); 5: Our work; 6:'Gao et al.' ('20011); 7: iGlenn & Hunter! ( l200lb ; 



: used for constra ining the models and for computing the C and CO cooling rates, 



Pour Overlap : 1 : iGerin & Phillips I lOOdti but spectra have been analy zed again; 2: Stanford et alj (Il990h ; 3: 
lAaho et al.' il99^; 4: 'Zhu et al.' (2003); 5: Our work; 6: 'Gao et alJ ('200ll); 7: iGlenn & Hunter! ( l200lb : *: used for 



constraining the models and for computi ng t he C and CO cooling rat es- 



Pour HENIZE 2-10 : 1: iBaas et alJ (Il99i): 2: !Kobulnickv et al] (!l995!) : 3: ?; 4: !Gerin & Phiuip3 (!2000!) but spectra 



have been analyzed again; 5: Bavet et all ( 2004!) 
CO cooUng rate s (seeiBavet et am2004 ). 
Pour ICIO : 1: 



used for constraining the models and for computing the C and 



Gerin & Phillips! ( |2000!) but spectra have been analyz ed a gain; 2: iBolatto et alj (l2000l) ; 3: Our 



work; 4: PhD Thesis of Becker (1990); 5: 'Petitpas & Wilson ('l998a|); 6: fclenn & Hunter! ( i200lh ; *: used for 
constraining the models and for comp uting the C and C O cooUng rate s.. 

Pour IC 342 : 1: llsrael & Baa3 (12003!) : 2: Our work; 3: lEckart et all (Il990l) : *: used for constraining the models 



and for computing the C and CO cooling rates. 



1: Petitpas & Wilson a998bi) : 2: llsrael & Baal ( l200ll) : 3: iLundgren et"ai] ( |2004 : 4: iHanda et al. 



Pour M 83 

(Il990l) : 5: Our wor k; *: used for constr ain ing the models an d fo r computing the C an d C O cooling rates , 



Pour NGC 253 : 1: 'Harris et al.' ('l99l"); 2: 'Wall et al.' ('l99l'); 



'Guesten et al.' (1993); 4: 



'Harrison et al.Vl995'); 



llsrael et al "(1995,): 6:.Mauersberger et al.(1996): 7: Harrison et al. (1999): 8: Dumke et al..(2001,) : 9: Israel & Baa3 
1 2002 ): 10: Bradford et al. 1 2003!): 11: iBavet et alj ( 12004 ): *: used for constraining the models and for computing 
the C and CO cooling rates (see'Bayet et al.' ('2004')). 

Pour NGC 6946 : 1: Israel & Baas, UOOl.) : 2: Gerin & PhimpsI (|2000!) but spectra have b een analyzed again; 3 : 
Walsh et alJ (12002!): 4:llshizuki et alJ 119901) : 5 :lweiiachew et alJ(ll988|): 6:ISofue et alJ(fT98i): 7:ICasoU et alJ(ll990!): 



8: Our work; 9: Iciausset et al.' (Il99l|); 10: |WaU et al.' (l99^; 11: |Mauersberger et alj (ll999|); 12: iNietoLeLaL 
(11999ft : *: used for constraining the models and for computing the C and CO cooUng rates. 



